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SUMMARY 

Tbe  equations  of  motion  of  an  aircraft,  for  small  perturbations  from  flight 
with  constant  linear  and  zero  angular  velocities,  are  developed  in  detail  using:- 

constant-velocity  or  body-fixed  axes, 
encastrd  or  free-free  modes, 

displacement  or  velocity  body  freedom  coordinates. 

The  relationship  is  clearly  stated  between  these  various  forms;  and  with  other 
proposed  forms,  in  particular  those  using  mean-body  axes.  The  whole  development 
is  kept  consistent,  as  far  as  possible,  with  the  Hopkin  notation  scheme 
(RAM  3562) . /'TvFhRTC  I 
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1 INTRODUCTION 

An  essential  feature  of  effective  research  is  good  communication.  For 
communication  to  be  successful  it  has  to  be  unambiguous,  and  both  attractive  and 
understandable  to  the  'recipient'.  One  wonders  how  much  good  work  has  been 
largely  wasted  not  because  it  has  not  been  reported  clearly  but  because  those 
trho  could  have  made  good  use  of  it  have  been  deterred  by  approaches  which  were 
strange  to  them  and  which  they  could  not  easily  relate  to  their  own  experience. 

Several  different  groups  of  aeronautical  research  scientists  have  been  led 
to  investigate  the  dynamics  of  aircraft  taking  some  account  of  their  flexibility. 
These  groups  include  those  whose  primary  concerns  are  active  control,  stability, 
ride  control,  structural  loads,  flutter  etc.  In  this  particular  field  ~ the 
dynamics  of  deformable  aircraft  - it  has  long  been  recognised  that  there  are 
serious  conmunication  barriers.  In  a recent  paper*  Taylor  and  Woodcock  each 
independently  sought  to  provide  a clear  statement  of  the  fundamentals  of  this 
subject.  The  present  paper  is  intended  as  a sequel  to  that  paper*  and  in 
particular  to  Part  II  of  that  paper*. 

The  objective  of  this  paper  is  to  study  in  detail  various  forms  of  the 
equations  of  motion,  for  small  perturbations  from  a datum  motion,  establish  the 
relationship  between  them,  and,  if  possible,  demonstrate  how  one  can  transform 
from  one  form  to  another. 

To  study  the  dynamics  of  deformable  aircraft  we  need  the  concepts  of  datum 
motion  of  the  aircraft,  and  the  undeformed  state  of  the  aircraft,  where  through- 
out the  datum  motion  the  aircraft  is  in  the  undeformed  state.  Departures  from 
the  datum  motion  are  called  perturbations.  Departures  from  the  undeformed  state 
are  called  deformations.  Consequently  we  need  two  fr2uiies  of  reference  - a 
datum-motion  frame  of  reference  and  an  undeformed-state  frame  of  reference  - 
which  are  such  that,  if  there  are  no  perturbations  the  aircraft  is  at  rest  rela- 
tive to  the  datum-motion  frame  of  reference,  and,  if  there  are  no  deformations 
the  aircraft  is  at  rest  relative  to  the  undeformed-state  frame  of  reference*. 

The  perturbations  comprise  therefore  (a)  the  translations  and  rotations  which  are 
necessary  to  make  the  datum-motion  frame  of  reference  coincide  with  the  undeformed- 
state  frame  of  reference,  and  (b)  the  deformations. 


* One  could  of  course  treat  all  perturbations  in  the  same  way  as  the  deformations, 
and  use  only  one  frame  of  reference  in  their  description,  but  from  several 
points  of  view,  and  in  particular  consideration  of  large  perturbations,  the 
above  procedure  seems  to  be  the  most  attractive. 
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Ue  restrict  ourselves  to  Lagrangian  methods.  Of  these  there  seen  to  be 
three  possible  types  of  approach: 


f 

t 

! 


(i)  The  use  of  Lagrange's  equation  for  an  inertial  frame  which  is  taken  as 
the  datum-motion  frame  of  reference. 

(ii)  The  use  of  Lagrange's  equation  for  a non-inertial  frame  (the  undeformed- 
state  frame  of  reference)  in  conjunction  with  the  principles  of  momentum 
equations  for  the  whole  aircraft. 

(iii)  The  same  as  (ii)  followed  by  some  coordinate  transformation. 

With  the  first  approach  there  is  usually  little  doubt  as  to  which  is  the  most 
appropriate  choice  of  axes  within  the  reference  frame.  The  axes,  within  that 
frame,  which  coincide  with  a particularly  significant  set  of  body-fixed  axes, 
such  as  the  principal  axes  of  inertia,  during  the  datum  motion,  would  be  a 
suitable  choice. 

When  the  second  approach  is  used  there  is  something  to  be  said  for  taking 
a non-inertial  frame  whose  position  etc  is  of  interest  in  itself.  We  have  there- 
fore taken  a certain  set  of  body-fixed  axes  to  define  the  non-inertial  frame. 

To  do  so  involves  imposing  certain  minor  restrictions  on  the  modes  of  deforma- 
tion. The  third  approach  can  be  thought  of  as  a way  of  removing,  to  all  intents 
and  purposes,  these  restrictions  at  the  expense  of  a loss  in  the  significance  of 
the  body-freedom  generalised  coordinates. 

The  notation  and  nomenclature  used  is  consistent,  apart  from  one  or  two 
noted  exceptions,  with  that  of  Fart  II  of  Ref  1,  and  therefore  also,  to  the  same 
degree, with  that  of  Hopkin's  comprehensive  scheme  . Attention  is  drawn  to  the 
Glossary  of  terms  and  List  of  symbols  at  the  end  of  this  paper.  We  confine  our- 
selves to  dimensional  forms  of  the  equations  of  motion  and  leave  any  normalisa- 
tion or  non-dimensionalisation  to  the  reader.  A further  paper*^,  which  is 
essentially  a particularisation  of  the  present  work  for  the  case  of  symnetric 
perturbations  in  heave,  pitch,  fore  and  aft  translation,  and  one  deformation  mode, 
does  however  also  give  suggested  non-dimensional  forms  of  the  equations. 

2 DATUM  MOTICW 

We  postulate  a datum  motion,  with  respect  to  which  the  dynamics  of  a 
deformable  aircraft  are  to  be  investigated.  One  particular  type  of  datum  motion 
is  considered:  straight  flight,  not  necessarily  level,  with  the  aircraft  having 
constant  linear  velocity  and  zero  angular  velocity.  It  is  assumed,  as  desired, 
that  no  deformation  of  the  aircraft  takes  place  during  the  datum  motion.  This 
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implies  that  all  the  forces  acting  on  the  aircraft  are  constant*  and  in  particu- 
lar, therefore,  that  the  atmosphere  is  uniform.  The  only  forces  acting  on  the 
aircraft  are  assumed  to  be  the  aerodynamic  and  propulsive,  gravitational  and 
structural  forces;  with  the  addition,  in  the  case  of  ground  contact,  of  what  we 
call  upholding  (or  support)  forces. 

As  a basic  frame  of  reference  a constant-velocity  axes  system  is  used  with 

origin  at  the  aircraft  centre  of  gravity  in  the  datum  motion  and  axes  parallel  to 

its  principal  axes  of  inertia  in  the  datum  motion.  They  are  therefore  body-fixed 

2 

axes  during  the  datum  motion;  in  fact  the  axes  which  Hopkin  calls  datum-attitude 
earth  axes.  Throughout  this  paper  we  will  use  the  term  constant-velocity  axes 
for  these  particular  axes.  The  subscript  f is  used  to  denote  values  of 
quantities  during  the  datum  motion. 

3 DEGREES  OF  FREEDOM 

The  aircraft  is  assumed  to  be  semi-rigid,  ie  having  a finite  number  of 
degrees  of  freedom  in  addition  to  its  six  body  freedoms.  Tr.o  representations  of 
any  possible  perturbations  are  used.  They  are  not  exactly  equivalent  but  to 
first  order  they  agree. 

The  first,  chosen  for  convenience  in  derivations  using  the  constant- 
velocity  axes  (alias  the  datum-attitude  earth  axes) , can  be  visualised  as 
follows.  The  transformation  of  the  aircraft  from  its  datum  motion  position  at 
any  instant  to  its  perturbed  position  and  shape  at  the  same  instant  can  be 
achieved  by  Che  following  successive  steps: 

(i)  Translations,  as  a rigid  body,  Xj  , y,  , - 
of  the  respective  constant-velocity  axes. 


1 


I 


in  the  directions 


(ii)  Successive  rotations,  as  a rigid  body,  i|),  0,  <ti  about  the  carried 
axes**  Oz,  Oy,  Ox  idiere  Oxyz  are  body-fixed  axes,  with  origin  at  the  particle 
(reference  point)  which  is  coincident  with  the  aircraft  centre  of  gravity  during 
the  datum  motion,  whose  orientation  is  fixed  in  a small  portion  of  the  aircraft 
which  includes  the  reference  point,  which  is  either  essentially  rigid  or  other- 
wise such  that  the  axes  will  always  remain  mutually  perpendicular.  The  orienta- 
tion of  the  axes  is  chosen  so  that  they  coincide  with  the  constant-velocity 
axes  during  the  datum  motion.  A frame  of  reference  which  coincides  with  the 

* We  are  assuming  also  that  the  aircraft's  mass  and  mass  distribution  are  con- 
stant and  so  neglecting  the  changes  produced  by  fuel  consumption. 

**  By  carried  axes  is  meant  the  position  of  the  axes  following  the  previous  rota- 
tions. It  is  important  to  remember  that  the  order  of  these  rotations  is  not 
commutative  in  general. 
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body-fixed  axes  after  this  step  will  be  called  the  no-deformation-body-fixed  axes. 
The  bracketed  superscript  n is  used  to  denote  reference  to  them.  They  are,  as 
will  be  seen,  in  general  only  body-fixed  as  long  as  the  aircraft  is  undeformed 
from  its  datum  shape.  It  should  be  noted  that  their  relevance  is  limited  to  this 
first  representation  of  the  perturbations. 

(iii)  Deformation  such  that  the  position  of  a particle  relative  to  the 
origin  of  the  no-deformation-body-fixed  axes,  and  referred  to  those  axes  is  given 
by 

r I 1 I 1 

(1) 


r (n)1 

X 

n 

..  * 

— 

*f 

+ R 

‘Jl 

• 

^n 

yf 

• 

• 

z 

n 

where  there  are  n deformational  degrees  of  freedom  represented  here  by  the  modal 
matrix  R whose  elements  are  functions  of  the  particle  being  considered  and  are 
independent  of  the  datum  motion.  For  a given  datum  motion  R will  be  a function 
of  (x^,yj,Zj).  The  function  R is  also  constrained  to  be  such  that  for  small 
perturbations  the  body-fixed  axes  remain  mutually  at  right  angles.  Writing 


R 


r 

..  b 


(2) 


I n 

. • • • c 

1 n 


the  direction  cosines  of  the  body-fixed  axes,  in  the  0 x y z (no-deformation- 

n n-'n  n 

body-fixed  axes)  reference  frame,  are  proportional  to 


and 


* O ’i  ■ ■ 

0 0 
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respectively;  and  so  they  remain  mutually  at  right  angles  for  small  perturbations 
provided  that 


©).  • ■©. 

ft).  - -ra. 

(S;).  • - 


for  all  i 


A more  complicated  deformational  representation  than  that  of  equation  (1),  or 
further  constraints  on  the  modal  matrix  R in  addition  to  (3),  would  be 
necessary  to  ensure  that  the  body^fixed  axes  remained  mutually  at  right  angles 
for  larger  perturbations. 

Thus,  with  this  first  representation  of  the  perturbation,  the  position  of 
a particle  relative  to  the  origin  of  the  constant-velocity  axes  and  referred  to 
those  axes  is  {of.  Appendix  A of  Part  II  of  Ref  1) 


where  S , a function  of  the  angles  6,  i{»  , is  the  axes  transformation  matrix 
(or  attitude  deviation  matrix)  of  Refs  I and  2.  Since 

S * I - A,  (5) 


where* 


r 0 -il)  el 


* Following  Refs  I and  2 the  symbol  A is  used  to  denote  a particular  skew- 
symnetric  matrix  formed  from  the  elements  of  a column  matrix  whose  leading 
element  is  used  as  a subscript  to  A . 


equation  (4)  can  be  written,  for  small  perturbations 


xJ  + IR  I -A^  ] q, 


where  q^^,. 


. , are  body  freedom  coordinates 


♦ . 


However  a better  approximation,  (54),  than  (7),  is  required  for  certain  purposes 
in  the  derivation  of  the  equations  of  motion  for  small  perturbations. 

The  alternative  representation  of  any  perturbations  from  the  datum  motion 
condition  is  to  define  the  deformation  relative  to  the  body-fixed  axes  (see 
Glossary  of  terms)  such  that  the  position  of  a particle  relative  to  the  reference 
point  is  given  by 

PI  ■ pfl  * ‘ "o  * 


where  is  the  value  of  R at  the  reference  point,  and  {of.  equation  (2)) 


Oc,/3yf)o  Oc^ayf)o 


(9.,/3zPo 


Ob./3x^)„ 


]] 


Remembering  that  the  elements  of  R satisfy  the  condition  (3),  it  is  easily  seen 
that  equation  (10)  always  represents  a perturbation  which  is  precisely  a deforma- 
tion relative  to  the  body-fixed  axes  - the  displacement  of  the  body-fixed  axes  is 
zero.  The  position  of  a particle  relative  to  the  origin  of  the  constant- 
velocity  axes  and  referred  to  those  axes  is  therefore  given  by 


(12) 


x<‘^> 

c 

- 

-(c) 

*1 

+ 
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+ (R  - 

-q,- 
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y(c) 

yf 
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-^f- 

in- 

> 

where  (x|*^\y are  the  translations,  and  are  the  rotations 

(according  to  the  standard  Euler  procedure)  which  transform  the  constant-velocity 
axes  into  the  body-fixed  axes.  Thus  for  small  perturbations  this  second  represen- 
tation gives 


(13) 
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It  will  be  seen  therefore  that  the  two  representations  agree  for  small  perturba- 
tions (equations  (7)  and  (13))  when  the  relationship  between  the  two  sets  of 
generalised  coordinates  is 


12 
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0 0 
I 0 
0 I 


'n+6 
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I 

0 

0 

-1 

’^1 

^0 

I 

0 

• 

• 
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p 

_ q 

0 

I 

*^n+6 

(16) 


There  is  however,  no  relationship  between  the  two  sets  of  generalised  coordinates 
which  will  ensure  agreement  of  the  two  representations  when  the  perturbations  are 
not  small  (ie  agreement  of  (4)  and  (12)). 

4 EQUATIONS  OF  EQUILIBRIUM 

The  constant-velocity  axes  (datum-attitude  earth  axes)  are  an  inertial 
frame,  since  the  datum  motion  has  zero  angular  velocity.  Consequently  the 
equations  of  motion  can  be  derived  from  the  inertial  frame  form  of  Lagrange's 
equation:— 


where  W is  the  kinetic  energy  relative  to  the  frame  of  reference,  and  the  Q^ 
are  the  total  generalised*  forces  in  the  degrees  of  freedom  obtained  by  the 
principle  of  virtual  work.  With  the  freedom  specified  by  equations  (4)  and  (7) 
it  is  easily  seen  that,  for  small  perturbations,  W is  a quadratic  form  in  the 
time  derivatives  q^  of  the  generalised  coordinates.  Consequently  the  existence 
of  the  datum  motion  requires  that  each  of  the  total  generalised  forces*  satisfies 
the  condition 

(Vf  = 0 (18) 

as  regards  their  unperturbed  values.  Now  from  (7) 


* We  have  referred  to  the  'total  generalised  force'  since  subsequently  we  describe 
the  various  constituents  as  'structural  generalised  force'  etc. 
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and  so  the  virtual  work  done  by  a distribution  of  force  vectors  (e^.f^.g^) 
referred  to  the  constant-velocity  axes,  acting  on  the  particle  whose  location  is 
(Xf,yf,Zf),  is,  in  the  datum  state. 


Thus 
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where  (X^.Y^.Z^),  (L^.M^.N^)  are  respectively  the  overall  forces  and  moments 
produced  by  the  above  mentioned  distribution  of  force  vecto:;i.  In  flight  the 
force  vector  on  any  particle  may  include  aerodynamic,  gravitational,  propulsive 
and  structural  constituents.  Other  constituents  such  as  the  magnetic  forces  may 
be  present  but  they  are  considered  to  be  negligible  as  regards  our  present 
purpose.  Thus  the  equation  (18)  can  be  rewritten  as  the  matrix  equation 

[-  (Qi)f  + (Gj)j  + (Pj)f  + (E^)^J  - 0 (23) 

where  the  expression  in  the  [ ] is  the  element  in  the  ith  row  of  the  matrix  and 
the  individual  terms  are  respectively  the  aerodynamic,  gravitational,  propulsive, 
and  structural  contributions  to  ■(Q£)£  • Expression  for  these,  derived  follow- 
ing (22)  and  with  the  assumptions  of  section  3.1  are  given  in  Table  1.  The 
structural  generalised  forces  in  particular  merit  some  comment.  The  structure 
cannot  exert  any  overall  force  or  moment  on  itself  - this  is  intuitively  true  and 
moreover  is  demonstrated  to  be  so  in  Appendix  A for  the  particular  case  of  an 
isotropic  elastic  body.  It  follows  therefore  that 


for  i “ (n+I) , . . . , (n+6) 


With  body-fixed  axes  we  can  as  demonstrated  in  Part  II  of  Ref  1 derive  the 
equations  of  motion  using  Lagrange's  equations  referred  to  a non-inertial  frame 
in  conjunction  with  equations  based  on  the  principle  of  momentum.  These  equations 
are  respectively 

where  Vq  is  the  centrifugal  potential  function,  W is  the  kinetic  energy 
relative  to  the  frame  of  reference,  G.  is  the  gyrostatic  force,  J.  is  a 
certain  coupling  force  between  rotational  body  freedoms  and  the  deformational 

A 

freedoms,  and  the  are  the  generalised  forces  obtained  regarding  the  frame  of 
reference  as  stationary  during  any  virtual  displacement  (of.  Ref  1,  Part  II);  and 
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where  X etc,  L etc,  are  respectively  the  resolutes  along  the  body-fixed  axes 

of  the  total  applied  force  and  torque  about  the  reference  point,  and  u etc, 

m 

are  the  resolutes  of  the  velocity  of  the  particle  at  (x,y,z).  For  perturbations 
from  our  assumed  datum  motion  the  angular  velocity  of  the  body-fixed  axes  is 
given  by  (see  Part  II,  Appendix  A of  Ref  1) 
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The  perturbed  position  of  a particle  will  be  taken  to  be  given  by  equation  (10) 
which  satisfies  the  condition,  necessary  for  the  application  of  equation  (25), 
that  the  position  of  the  reference  point  and  the  orientation  of  the  body-fixed 
axes  is  independent  of  the  degrees  of  freedom  q.  (i  > 1,  ...,  n). 


Since  the  linear  velocity  is  constant  and  the  angular  velocity  is  zero 
during  the  datum  motion  it  is  found  (see  section  6.3  of  Part  II  of  Ref  1)  that 
all  the  terms  on  the  left-hand  side  of  equation  (25)  and  on  the  right-hand  sides 
of  equations  (26)  and  (27)  are  zero  when  there  are  no  perturbations.  Conse- 
quently the  datum  motion  equilibrium  conditions  are 
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It  is  easily  shovm  (of.  equation  (10))  that 
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which  can  be  alternatively  written  as 
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We  can  as  before  express  these  generalised  forces  as  a sum  of  aerodynamic, 
propulsive,  gravitational  and  structural  contributions.  Expressions  for  these  are 
given  in  Table  2. 


Finally  we  note  that  the  second  form  of  the  datum  motion  equations  (equations 
(29),  (30)  and  (31))  is  merely  the  first  form  (equation  (18))  premultiplied  by 
the  matrix 
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0 
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(34) 


5 APPLIED  FORCES 

5. 1 Assumptions 

The  primary  object  of  this  study  is  the  aircraft  in  flight  well  away  from 
the  ground.  When  the  aircraft  is  flying  near  the  ground  certain  assumptions  in 
respect  of  the  form  of  the  aerodynamic  forces  cannot  be  made,  while  when  it  is 
actually  in  contact  with  the  ground  additional  upholding  (or  support)  forces  have 
to  be  taken  into  account.  These  two  cases  are  considered  in  Appendices  B and  C 
respectively  while  here  we  restrict  ourselves  to  situations  where  there  is  no 
effect  at  all  of  the  near  presence  of  the  earth  other  than  gravity  and  the  fact 
that  the  flight  is  atmospheric. 

5.2  Aerodynamic 

The  atmosphere  is  assumed  to  be  uniform,  and  the  aerodynamic  forces  are 
assumed  to  have  no  hereditary  constituent.  Both  of  these  are  obvious  departures 
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from  the  truth.  The  first  will  rarely  produce  much  error  provided  the  'uniform' 
atmosphere  has  the  same  properties  - density  etc  - as  the  actual  atmosphere  at 
some  mean  altitude  of  the  aircraft's  motion.  For  a perturbation  which  is  a 
maintained  oscillation  the  second  assumption  is  in  effect  eliminated  by  allowing 

the  aerodynamic  coefficients  to  be  functions  of  the  frequency  of  the  oscillation. 

4 5 

Moreover  Woodcock  and  Lawrence  have  shown  ’ that  a good  approximation  to  any 
calculated  oscillatory  motion  (decaying,  maintained  or  growing)  is  obtained  pro- 
vided the  aerodynamics  used  is  for  maintained  oscillations  of  the  same  frequency. 
In  such  cases  one  assumes  a frequency  for  the  aerodynamics,  calculates  the  motion 
and  then  repeats  the  procedure  until  the  assumed  and  calculated  frequencies  are 
in  sufficient  agreement. 

The  boundary  condition  to  be  satisfied  by  the  air  motion  at  the  surface  of 
the  aircraft  involves  the  velocity  of  the  surface  normal  to  itself.  With  the 
perturbation  representation  of  equation  (10)  the  slope  vector  of  the  body  surface, 
referred  to  the  body-fixed  axes,  is  a function  only  of  the  generalised 
coordinates  qj,  ...,  q^  . Also  the  velocity  of  a particle,  referred  to  the  body- 
fixed  axes,  is  (see  Ref  1,  Part  II,  section  6.2  and  Appendix  A) 


- (R  - Rq  * Wq^  q,  + s Uj  + ij* 
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where  the  angular  velocity  vector  of  the  body-fixed  axes  is 


Pi  ♦ 


and  {x,y,z}  is  given  by  equation  (10). 


Consequently,  having  assumed  no  ground  and  no  hereditary  effects,  the  local 
aerodynamic  force  vector,  referred  to  the  body-fixed  axes,  will  have  the  form, 
for  small  perturbations: 


(37) 


Similarly  we  find  that  with  the  alternative  representation  of  the  deforma- 
tion the  local  aerodynamic  force  vector,  referred  to  the  constant-velocity  axes, 
will  have  the  form,  for  small  perturbations: 

- J 


* Terms  such  as  downwash,  sidewash  are  commonly  used.  This  is  a generalisation 
to  suit  any  surface. 
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Note  that  the  coefficients  ej,  above  expressions  (37) 

and  (38)  are  in  general  differential  operators  (polynomials  in  the  differential 
operator  D ■ d/dt). 

5.3  Gravitational 

One  can,  for  all  intents  and  purposes*,  assume  the  gravitational  force 

acting  on  a particle  of  mass  5m  to  be  a constant  force  5mg  acting  in  the  z 

2 ^ 

direction  (vertically  downwards)  of  the  normal  earth-fixed  axes  . The  orienta- 
tion transformation  from  these  axes  to  the  constant-velocity  axes  can  be 
achieved  by  a standard  Euler  sequence  of  successive  rotations  about 

carried  axes.  Thus  if  we  denote  the  last  column  of  the  axes  transformation 
matrix  S by  the  symbol  , then  the  local  gravitational  force  vector 
referred  to  the  constant-velocity  axes  is  5mgl^^  where 

= - sin  e . (39) 

f ^ 

sin  cos  6^ 

cos  cos  6^ 

Referred  to  body-fixed  axes  the  vector  components  are  not  constant  but  have  some 
terms  which  depend  on  the  perturbations.  Thus  we  have 


* The  rate  of  change  of  g with  location  is  extremely  small  almost  everywhere. 


5.4  Propulsive 

A very  simple  model  of  the  propulsive  forces  has  been  assumed.  Time  no 
doubt  will  tell  whether  a more  sophisticated  one  is  required.  It  is  assumed 
that  the  propulsive  force  acting  on  any  particular  particle  has  constant  compo- 
nents »^pf »8p£)  in  the  direction  of  the  body-fixed  axes. 

Since  the  modal  matrix  R (equations  (1)  and  (2))  satisfies  the  condition 
(3) , we  have 
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where  P^  is  given  by  equation  (11);  and  so,  with  the  deformation  of  equation  (1) 
the  position  of  a particle  relative  to  the  reference  point  (the  origin  of  the 
body-fixed  axes)  and  referred  to  the  no-deformation-body-fixed  axes  is,  for 
(Xf,y£,Zj)  small. 
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(42) 


Thus,  for  small  perturbations,  the  orientation  transformation  from  the  body- 
fixed  axes  to  the  no-deformation-body-fixed  axes  is  achieved  by  the  rotations* 
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'"q, 


l5nj 


(43) 


The  last  term  in  (42)  is  the  first  approximation  to  (S  - Dpq  where  S is 


the  appropriate  transformation  matrix  (.of.  equation  (5)).  Consequently  the 
orientation  transformation  from  the  body-fixed  axes  to  the  constant-velocity 
axes  is  the  result,  for  small  perturbations  of  the  rotation 
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* The  order  is  iomaterial  since  the  rotations  are  small. 
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We  find  therefore  that  the  local  propulsive  force  vector,  referred  to  the 
constant-velocity  axes,  has,  for  the  deformation  of  equation  (4),  the  value 
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5.5  Structural 


On  a particle  of  the  aircraft  it  is  assumed  there  is  what  we  call  a 
structural  force  produced  by  the  stresses  in  the  material  and  possibly  also 
other  causes  such  as  friction.  This  structural  force  is  taken  to  be  entirely 
determined  by  the  current  shape  of  the  aircraft.  Thus  the  two  forms  of  the 
structural  force  vector,  referred  respectively  to  the  constant-velocity  axes  and 
to  the  body-fixed  axes,  that  we  will  use  are;— 
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The  various  coefficients  in  the  above  expression  are  taken  to  be  constant. 
That  is  we  are  assuming  there  is  no  structural  damping.  No  really  adequate 
theory  of  the  damping  of  aircraft  structures  has  been  developed.  One  either 
neglects  it  and  assumes  that  by  doing  so  one  will  be  'on  the  safe  side'  - and  this 
seems  to  be  at  least  nearly  always  so;  or  else  one  includes  coefficients  in  the 
final  equations  based  on  ground  resonance  test  measurements  and  primitive  theory. 
It  may  well  be  that  for  most,  if  not  all,  aircraft  dynamic  studies  one  or  other 
of  these  procedures  is  adequate.  It  will  be  noted  that  we  could  here  fairly 
easily  introduce  a simple  representation  of  structural  damping  by  replacing  the 
constant  coefficients  of  equations  (46)  and  (47)  by  differential  operators. 

What  exactly  is  meant  by  the  term  'structural  force'  is  perhaps  made  clearer 
by  the  analysis  of  Appendix  A for  an  isotropic  elastic  body.  Intuitively  one 
knows  that  the  structural  forces  on  the  aircraft  will  not  produce  any  overall 
force  or  moment  on  the  whole  aircraft  (of.  Appendix  A) , but  there  may  well  be 
forces  tending  to  deform  the  aircraft. 

It  should  perhaps  be  emphasised  that  the  structural  force  coefficients 
e . etc  of  equations  (46)  and  (47)  , in  addition  to  the  datum  motion  structural 

S X 

force,  will  in  general  depend  on  the  chosen  unperturbed  condition  (of  equation 
(A-13)).  Thus,  for  example,  the  coefficients  for  the  ground  resonance  test 
situation  will  be  different  from  those  for  the  aircraft  in  level  flight  at  a 
given  speed.  Nevertheless  one  will  often  have  to  assume,  particularly  when 
structural  force  information  is  obtained  from  experiments  on  the  ground 
(of.  Appendix  C),  that  the  structural  force  coefficients  are  independent  of  the 
chosen  unperturbed  condition.  In  such  circumstances  there  is  usually  little  hope 
of  doing  anything  better. 

6 DERIVATION  OF  THE  EQUATIONS  OF  MOTION  FOR  SMALL  PERTURBATIONS 
6. I Equations  of  motion  using  constant-velocity  axes 

With  our  chosen  constant-velocity  axes  the  equations  of  motion,  as  already 
pointed  out  in  section  4,  can  be  derived  using  the  inertial  form  of  Lagrange's 
equation: 

_d_  Aw  \ _ ^ 
dt  ~ 
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The  kinetic  energy  W , relative  to  the  frame  of  reference,  is,  using  (7)  given 
by 


W - iX6mt4,  ...  Va’  r 1 ^ -V  r ^ 


and  so 
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1 6m  r’^R  R^  -r\  q 

Xf  j 


A R A -(A  )‘  q , 

x^  x^  x^  I ljn+6 


Since,  during  the  datum  motion,  the  reference  point  is  the  centre  of  gravity 
and  the  constant-velocity  axes  are  the  principal  axes  of  inertia,  we  have 

[6mA  » 0 (52) 

*f 

and  ^ 

- [6m(Ax^)^  = diagCl^  1^  I^}  = (53) 

where  I , I , I are  the  principal  moments  of  inertia. 

X y Z 

To  determine  the  virtual  work,  and  hence  the  generalised  forces,  we  need 
to  take  a better  approximation,  to  the  expression  for  the  position  of  a particle, 
than  that  (equation  (7))  which  sufficed  for  the  unperturbed  state.  From  equation 
(4) 
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* We  have  B . * i (^a  * ^a*)  where  the  subscript  dip  denotes  the  A matrix 
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Consequently, 


077 


' I 

I f 


where  the  expressions  for  the  (Qj^)£  have  already  been  given  (equation  (22)) 
and  where 
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This  expression  can  be  simplified  for 
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and  (.of  equation  (57))* 
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where  is  the  lower  triangular  matrix,  with  zeroes  on  the  principal  diagonal, 
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It  is  easily  be  seen  that 


Certain  of  the  above  terms  can  be  expressed  in  terms  of  translational  and 
rotational  force  coefficients.  The  overall  translational  forces  on  the  aircraft 
referred  to  the  constant-velocity  axes,  produced  by  the  distribution  of  local 
force  vectors  is  given  by  (of  equation  (65)) 
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The  overall  moments  about  the  constant-velocity  axes  produced  by  the  same  force 
distribution  are*  (of  equations  (7)  and  (65)) 
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* A subscript  c is  required  with  the  symbol  for  the  overall  moments  but  not  for 
the  overall  forces  since  only  the  former  depends  on  the  position  as  well  as  the 
orientation  of  the  constant-velocity  axes. 
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At  this  point  it  is  worth  noting  that  the  perturbation  in  the  moments  is  not  just 
the  moments  produced  by  the  perturbations  in  the  local  forces,  that  is  not  just 
the  second  term  in  the  above  expression.  It  should  also  be  clearly  understood 
that  the  moment  we  are  considering  is  the  moment  about  the  origin  of  the  constant- 
velocity  axes  and  not  about  the  reference  point  - hence  the  use  of  the  subscript 
c in  the  symbols  for  the  resolutes  of  this  moment  (eg  (76)). 


Thus  finally  we  have 
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Equations  (50)  (51)  and  (78)  along  with  (48)  provide  the  means  for  writing 
out  in  full  the  equations  of  motion  for  small  perturbations  from  the  datum  motion. 
These  are  given  in  Table  3.  The  various  contributions,  aerodynamic  etc,  to  the 
generalised  coefficients,  have  been  obtained  using  the  local  load  vectors  given 
in  sections  5.2  to  5.5  (equations  (38),  (39),  (45)  and  (46)).  One  may  be 
surprised  that  the  expression  for  the  propulsive  matrix  [P..l  given  in  Table  3 
involves  the  constant  matrix  (see  equation  (11)).  This  is  a consequence  of 

our  assumption  that  the  propulsive  force  acting  on  any  particular  particle  has 
I constant  components  in  the  direction  of  the  body-fixed  axes. 

j 6.2  Equations  of  motion  using  body-fixed  axes 

Initially  we  will  derive  the  equations  of  motion,  using  body-fixed  axes, 
from  Lagranges  equation  for  a non-inertial  frame  (equation  (25))  and  the  equations 
based  on  the  principle  of  momentum  (equations  (26)  and  (27))  given  in  section  4. 
The  perturbed  position  of  a particle  relative  to  the  body-fixed  axes  is  that  given 
in  equation  (10),  and  for  the  body  freedom  coordinates  we  take  the  translations 
{x|*^\  and  the  rotations  ($,0,'$)  (c/  equations  (14)  and  (15)).  With 

this  representation  the  deformational  freedcmis  are  deformations  relative  to  the 
j body-fixed  axes,  which  we  will  call  encastre  modes.  Thus  the  initial  derivation 

I , 

' can  be  described  as  for: 

(i)  encastr6  modes  and  displacement  body  freedoms. 

In  the  equations  thus  derived  we  can  make  coordinate  transformations  which 
will  give  different  meanings  to  the  generalised  coordinates,  and  so  from  the 
initial  equations  we  derive  the  equations  of  motion  in  terms  of: 

(ii)  encastr6  modes  and  velocity  body  freedoms; 

(iii)  free-free  modes  and  displacement  body  freedoms;  and 

(iv)  free-free  modes  and  velocity  body  freedoms. 

6.2.1  Equations  of  motion  in  terms  of  encastr6  modes  and  displacement 
body  freedoms 

The  determination  of  what  we  call  the  ponderous  terms,  those  not  resulting 
from  any  force*  but  from  the  fact  that  a system  having  mass  is  moving,  has  been 
considered  in  some  detail  in  Ref  1 (Part  II) . It  is  there  shown  that  the 
ponderous  dampings  and  stiffnesses  will  all  be  zero  if  the  angular  velocity  of 
the  aircraft  is  zero  throughout  the  datum  motion.  Indeed,  for  our  chosen  datum 

' * One  does  of  course  speak  of  such  things  as  inertia  forces,  centrifugal  forces, 

Coriolis  forces  but  these  are  not  really  forces  but  merely  a convenient 
^ analogy  arising  from  D'Alembert's  principle. 

LL_j 


077 


33 


motion,  in  the  non-inertial  form  of  Lagrange's  equation  (see  Ref 
Ref  15), 


Part  II  or 


(79) 


there  are  no  first  order  terms,  in  the  generalised  coordinates,  arising  from 
or  3W/3q.  . 

The  linear  velocity  of  the  reference  point  and  the  angular  velocity  of  the 
body-fixed  axes  are  (of  Ref  1,  Part  II  section  6.4)  given  by 
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Thus,  with  a particle  position  relative  to  reference  point  defined  by  equation 
(10),  the  velocity  of  the  particle  is  given  by 
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To  the  same  order  of  approximation  the  velocities  of  the  reference  frame 
(equations  (80)  and  (81))  are 


The  centrifugal  potential  function  Vq  is  therefore  for  small  perturbations 
given  by 


and  so 


The  J.  terms  are  given  by 
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The  kinetic  energy  relative  to  the  frame  of  reference  (the  body-fixed  axes)  is 
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and  so 
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The  above  equations  (86) , (87)  and  (89)  suffice  to  give  the  ponderous  inertia 
coefficients  for  i - 1 ->■  n and  all  j (see  equations  (14)  and  (15)). 
Details  are  given  in  Table  4. 


The  equations  based  on  the  principle  momentum  are,  from  (26),  (27)  and  (82) 


and  so  the  expressions  on  the  right  of  these  immediately  give  the  other  ponderous 
inertia  coefficients  (of  Table  4) . 

Due  to  the  force  vector 


(92) 


it  is  easily  seen  that  the  translational  force  and  moment  about  the  body-fixed 
axes  are,  respectively 


(93) 


and 

y 


Details  of  the  separate  contributions  (aerodynamic  etc)  to  these  forces,  obtained 
using  the  different  local  force  vectors  (equations  (37) , (40) , (47)  and  section 
5.4),  are  given  in  Table  9.  As  pointed  out  in  section  5.5  and  Appendix  A,  there 
will  be  no  structural  contribution. 

For  the  deformational  freedoms  the  generalised  forces  resulting  from  the 
local  force  vectors  (92),  are,  using  (10)  and  remembering  that  the  frame  of 
reference  is  regarded  as  stationary  in  the  assumed  virtual  displacement  used  to 
calculate  the  virtual  work. 
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The  individual  contributions  (aerodynamics  etc)  to  these  generalised  forces, 
obtained  using  the  different  local  force  vectors  (equations  (37),  (40),  (47)  and 
section  5.4),  are  given  in  detail  in  Table  4. 

We  have  therefore  obtained  the  complete  set  of  equations  of  motion  for 
small  perturbations  of  the  datum  motion.  The  ponderous  terms  come  from 
equations  (86),  (87),  (89),  (90)  and  (91),  and  the  various  forces  from  equations 
(93),  (94)  and  (96).  Table  4 gives  the  equations  written  out  in  full. 

It  may  be  remarked  that  anybody  coming  to  use  the  expressions  of  Table  4 

may  well  not  know  the  individual  matrices  R,  R^,  but  only  the  complete  modal 

matrix  (R  - Rq  + A^^P^)  {of  equation  (10)).  The  constant  matrices  Rq  and  P^ 

are,  however,  in  a sense  arbitrary.  They  have  only  been  introduced  to  ease  the 

understanding  of  the  relationship  between  one  form  of  the  equations  of  motion 

and  another.  Thus  if  constant  terms  are  added  to  R„  and  P , and  correspond- 

0 q 

ing  changes  made  to  R so  that  the  new  matrix  (R  - R^  + Ax^P^)  is  exactly  the 
same  as  the  original  one,  the  new  matrix  R will  still  satisfy  the  condition 
(3).  To  set  up  the  equations  of  motion  using  Table  4 therefore  one  might  as 
well  take  P^  and  Rq  to  be  zero,  but  it  must  be  remembered  that  R must  then, 
of  course,  be  a matrix  of  encastr€  modes. 

6.2.2  Equations  of  motion  in  terms  of  encastre  modes  and  velocity  body 
freedoms 

The  linear  velocity  of  the  reference  point  and  the  angular  velocity  of  the 
body-fixed  axes  are  given  by  equations  (80)  and  (81).  Taking  the  perturbation  of 
these  components  from  their  datum  motion  value  as  new  generalised  coordinates  for 
the  body  freedoms  we  have 
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and 
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The  generalised  coordinates  for  the  deformational  freedoms  are  kept  the  same  as 
those  used  in  section  6.2.1  and  so  all  together  we  have  the  coordinate 
transformation 


■ 

• 

3n+6_ 

0 ID  Au 


f 

ID 


‘1 


L?n+6J 


(99) 


which  can  be  inverted  to  give 
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(100) 


since  both  sets  of  coordinates  will  be  zero  in  the  unperturbed  state. 

One  could  start  afresh  and  obtain  the  equations  of  motion,  as  in  section 
6.2.1,  using  the  new  coordinates.  It  is  obvious  however  that  the  same  result  is 
achieved  if  the  transformation  (100)  is  applied  to  the  equation  of  motion  for  the 
case  (encastrd  modes,  displacement  body  freedoms)  of  section  6.2.1  and  Table  4. 
This  has  been  done  with  the  outcome  shown  (however  see  below)  in  detail  in 
Table  5.  It  will  be  seen  that  the  ponderous  inertias  now  give  rise  to  ponderous 
dampings  and  ponderous  stiffnesses  . One  would  similarly  expect 

terms  in  D"!  and  D ^ from  the  various  generalised  force  matrices.  However 
because  of  the  form  of  these  matrices,  and  in  particular  the  form  of  the  aero- 
dynamic matrix  for  flight  well  away  from  the  ground,  the  only  negative  power  of 
the  differential  operator  D is  that  associated  with  the  gravitational  force 
(see  Table  5) . ^ 

It  may  have  been  considered  desirable  to  have  3 symmetric  inertia  matrix. 
However  this  is  not  possible  - to  make  it  symmetric  one  finds  that  one  has  to 
premultiply  the  equations  of  motion  by  a singular  matrix  and  so  would  just  throw 


40 


away  some  of  our  information.  The  submatrix  which  is  the  inertia  matrix  for  the 
degrees  of  freedom  1 ->•  n (the  deformational  freedom)  will  however  be  symnetric. 

In  Table  5 some  new  symbols  are  used  resulting  from  writing  the  local  aero- 
dynamic force  vector  (equation  (37))  in  terms  of  the  new  generalised  coordinates 
(equation  (99)),  ie 


Z pj]  ’3  * E pj] 
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(101) 


where 


i = 1 ...  6 . (102) 


In  the  use  of  Table  5 one  can,  as  pointed  out  at  the  end  of  section  6.2.1, 
take  the  matrices  P and  R„  to  be  zero  and  R to  be  a matrix  of  encastre 

q 0 

modes.  The  form  given  there  is  consequent  upon  the  consistent  use  of  R through- 
out the  paper  for  a matrix  of  free-free  modes. 

6.2.3  Equations  of  motion  in  terms  of  free-free  modes  and  displacement 
body  freedoms 

As  already  pointed  out,  the  use  of  the  non-inertial  form  of  Lagrange's 
equation  (equation  (79))  necessitates  the  choice  of  coordinates  for  the  deforma- 
tional degrees  of  freedom  such  that  the  position  and  orientation  of  the  reference 
frame  (the  body-fixed  axes  in  our  case)  are  independent  of  these  coordinates. 
However  one  may  often  prefer  to  have  deformational  freedoms  which  involve  some 
movement  of  the  body-fixed  axes,  as,  for  example,  when  one  desires  to  use  normal 
modes  as  some  of  the  degrees  of  freedom.  This  difficulty  is  simply  overcome  by 
applying  a coordinate  transformation  to  the  equation  of  motion  obtained  in 
section  6.2.1. 


If  we  put 
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(104) 


then  the  total  perturbation  can  be  thought  of  as  made  up  of  the  successive 
perturbations : 

(i)  Translations,  in  the  direction  of  the  constant-velocity  axes,  as  a rigid 
body: 
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(ii)  Rotations  as  a rigid  body,  according  to  the  standard  Euler  procedure: 
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(iii)  Deformations: 
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The  position  of  a particle  relative  to  the  origin  of  the  constant- velocity  axes 
and  referred  to  those  axes  is  therefore,  from  (12),  for  small  perturbations: 
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Each  of  the  generalised  coordinates  qj  "*■  q is  a measure  of  perturbation  just 
in  an  unconstrained,  or  free-free  mode  - there  are  no  restrictions*  on  R as 
distinct  from  (R  - Rq  ♦ AxfPq) . 

As  well  as  making  the  coordinate  transformation  (103)  we  premultiply  the 
equation  of  motion  by  the  matrix 
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so  that  the  symmetry  of  the  ponderous  inertia  matrix  will  be  preserved.  The 
resultant  equation  of  motion  is  given  in  detail  in  Table  6. 

The  local  aerodynamic  force  vector,  in  this  case,  is  written  as  (.of 
equation  (37)) 
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where  the  aircraft  is  well  away  from  the  ground. 

6.2.4  Equations  of  motion  in  terms  of  free-free  modes  and  velocity  body 
freedoms 

To  change  from  displacement  to  velocity  generalised  coordinates  we  made  the 
transformation  (99) . To  change  from  encastr^  modes  to  free— free  modes  we  made 
the  transformation  (103).  To  make  both  these  changes  we  do  the  latter  followed 
by  the  former  and  so  we  put 


* Apart  of  course  from  the  general  restriction  that  we  have  imposed  throughout 
that  it  is  such  that  the  body-fixed  axes  always  remain  an  orthogonal  frame. 
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The  linear  and  angular  velocities  of  the  body-fixed  axes  are  therefore,  from 
(97)  and  (98) 
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Thus  the  coordinates  "*■  ^^+6  components  of  the  linear  and  angular 

velocities  of  the  body-fixed  axes  relative  to  their  values  due  to  the  deformation 
in  the  free-free  modes. 

The  position  of  a particle  relative  to  the  origin  of  the  constant-velocity 
axes  and  referred  to  those  axes  is  from  (12),  for  small  perturbations: 
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whereas  with  the  encastr^  modes  of  section  6.2.2  it  was,  from  equations  (12)  and  ^ 

(100) , it  was :- 
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In  the  present  case  therefore  the  generalised  coordinates  5]  measure 

perturbations  in  unconstrained  (free-free)  modes,  in  contrast  to  the  'encastrd 
modes'  which  were  constrained  to  have  zero  displacement  and  slope  at  the 
reference  point. 

It  is  desirable  that,  for  the  deformational  freedoms  alone,  if  not  for  all 
the  freedoms,  the  inertia  matrix  should  be  symmetric.  This  is  so  for  all  the 
previous  developments  and  we  arrange  it  to  be  so  now  by  premultiplying  the 
equations  of  motion  by 
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in  addition  to  the  coordinate  transformation  (112).  The  resultant  equation  of 
motion  is  given  in  detail  in  Table  7. 

In  terms  of  the  current  coordinates  the  local  aerodynamic  force  vector  is 
written  as  (see  equation  (37)) 
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1 ...  6 


(119) 


the  aircraft  oeing  well  away  from  the  ground. 
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(121) 


* The  dressings  refer  respectively  to  the  cases  encastre  modes- 

displacement  body  freedoms,  encastrd-velocity,  free-f ree-displacement,  and 
f ree-f ree-ve loc i ty . 
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Referring  the  latter  to  the  constant-velocity  (ie  datum-attitude  earth)  axes 
we  have 
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where,  from  equations  (98),  (104)  and  (112) 
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If  the  typical  particle  position  vectors,  relative  to  the  origin  of  the  constant- 
velocity  axes  and  referred  to  the  same  axes  are,  to  first  order  of  smallness,  the 
same  in  each  derivation  then  the  force  vectors  will  be  the  same,  to  first  order, 
in  each  derivation.  Thus  if  (of  equations  (7)  and  (13)) 
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then  the  corresponding  force  vectors  will  be  equal  for  small  perturbations. 
Equation  (124)  is  satisfied  if 
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For  the  other  sets  of  generalised  coordinates  (124)  is  satisfied  by  (see 
equations  (99),  (103)  and  (111)) 


Equations  (126)  and  (130)  are  the  required  relationships  between  the  various 

forms  of  the  local  force  vector.  Putting  them  another  way  they  are 
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7.1.1  Between  local  aerodynamic  force  vectors  when  well  away  from  the 
ground 


For  an  aircraft  in  flight  well  away  from  the  ground,  as  in  the  assuiiq>tion 
throughout  the  main  part  of  this  paper,  it  can  reasonably  be  assumed  that  an 
adequate  approximation  to  the  local  aerodynamic  force  vector  has  the  form  given 
in  equation  (37)  (for  body-fixed  axes,  encastrd  modes  and  displacement  body 
freedoms)  or  (38)  (for  constant-velocity  axes).  That  is,  with  expressions  of 
the  forms  given  in  equations  (120)  or  (121),  certain  coefficients  have  the 
following  character: 
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It  follows  from  (126)  that 


With  the  other  sets  of  generalised  coordinates  we  have  (af  equations  (101),  (102) 
(110),  (118)  and  (119)) 
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It  will  be  noticed  that  the  relationships  (139),  (141)  and  (142)  have  already 
been  substituted  in  equations  (101)  and  (110),  as  appropriate,  and  the  symbols 


e . 
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(i  =■  1 ->■  n)  and  e"  . . (i  = I 6)  are  not  used  outside  this  section. 
n+1 


Equation  (138)  gives  the  relationship  one  would  expect  between  the  two  cases 
where  one  has  free-free  deformation  modes  (the  ~ and  the  *');  (139)  what  is 
expected  for  the  two  cases  with  encastr^  modes  (the  ''  and  the  (140)  what  is 

expected  for  the  two  cases  with  velocity  rigid  body  coordinates  (the  ''  and 
the  '' ) ; and  (141)  what  is  expected  for  the  cases  with  displacement  rigid  body 
coordinates  (the  * and  the  ~) . Finally  it  is  interesting  to  note  that,  from 
(136)  and  (138), 
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This  relationship  could,  of  course,  have  been  immediately  obtained  from  (130). 
Equations  (139)  and  (143)  are  in  fact  independent  of  the  particular  form  of  the 
aerodynamic  forces  consequent  upon  remoteness  from  the  ground. 
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7.2 


Between  overall  force  and  moment  vector  coefficients 


From  the  expressions  for  the  overall  forces  and  moments  (equations  (74), 
(76),  (93)  and  (94))  and  the  relationship  (131)  between  the  local  force  vectors 
it  is  easily  seen  that  we  have  the  following  connection  between  the  coefficients 
for  the  constant-velocity  axes  development  and  those  for  the  basic  bodjr-fixed 
axes  (encastre  modes  displacement  body  freedoms)  development: 
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Using  the  relationships  between  the  various  sets  of  body-fixed  generalised 
coordinates  (equations  (99),  (103)  and  (111)),  we  also  immediately  find  that: 
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with  the  particular  forms  of  the  local  aerodynamic  force  vectors 
appropriate  to  flight  where  there  is  no  ground  effect  (equations  (37) , (38) , 
(101),  (110)  and  (118))  it  follows  that  certain  coefficients  have  particular 
forms  {of  Tables  8 to  12) 


and  similarly  for  the  ^ and  ~ coefficients  except  that  the  term  Ax^,  Al^ 
are  absent.  For  the  **  coefficients  (and  similarly  for  the  ''  coefficients) 
we  can  simply  write  X , * X , X , * X etc.  It  is  then  easily  seen 


(of  section  7.2)  that 
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The  above  relationships  can  be  confirmed  from  the  expression  of  Tables  8 to  12 
using  the  relationships  between  the  local  aerodynamic  force  vectors  given  in 
section  4.1.1.  Anybody  with  prior  knowledge  of  relationships  between  aero- 
dynamic derivatives  may  well  be  surprised  by  the  equations  (152).  They  are  used 
to  relationships  such  as  (see  eg  Ref  2,  section  22),  for  our  case  of  zero 
angular  velocity  during  the  datum  motion. 


Z + w,Z*  - v^Z* 
p f V f w 


(157) 


It  must  be  emphasised  therefore  that  our  terms  Z*  are  not  derivatives.  In 

<P 

general  they  are  differential  operators.  Thus  we  could  write  (the  notation  is 
ephemeral) 
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and  then,  making  use  of  (97)  and  (98),  we  find  that  the  coefficients  of 

A M 

Z and  Z are  respectively 


in 


and 


Z:  + Z..W-  - z„v. 
Ip  y f 2 f 


Z + Z*W,  - Z'V^ 
p V f w f 


(159) 


The  former  term  is  commonly  called  Z^  since  it  has  come  from  a Taylor  expansion 

in  terms  of  the  displacements;  while  the  latter  term  is  written  as  the  linear 

( 

combination  (Z^  + Z^w^  - Z^v^)  of  three  derivatives  since  a Taylor  expansion  in 
terms  of  velocities  was  first  used  before  transforming  to  the  displacement 
coordinates.  Equating  these  two  expressions  then  gives  (157)  which  contrasts 
with  our  relationship 

Z-  - Zp  (160) 

where  the  quantities  involved  are  not  derivatives. 

7 • 3 Between  coefficients  in  perturbation  equations  of  motion 

In  the  constant-velocity  axes  derivation  a general  expression  [Q. .]  for 

any  of  the  matrices,  [Q.  , -[E,.],  -tPii],  -[G..]  was  obtained  in  equation 

ij  ij 

(78).  The  corresponding  general  matrix  [Q^^]  for  the  body-fixed  axes  derivation 
using  encastre  modes  and  displacement  body  freedoms  is  given  by  equations  (93) , 
(94)  and  (96) . From  the  relationships  between  the  local  force  vectors 
(equation  (131))  and  between  the  overall  forces  (equations  (144)  and  (145))  we 
then  find  that 


I 


In  addition  we  have  the  relationship  between  the  ponderous  inertia  matrices 
(of  Tables  3 and  4) 


The  three  other  forms  of  the  perturbation  equation  of  motion  that  we  have  con- 
sidered were  all  obtained  simply  from  the  equation  for  the  body-fixed  axes, 
encastr4  modes,  displacement  body  freedoms,  case  by  a coordinate  transformation 
and  possibly  also  a matrix  premultiplication  (see  sections  6.2.2,  6.2.3  and 
6.2.4).  Thus  we  have 
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where  each  of  these  matrices  may  be  a Q. . , -E. . , -P. . , or  -G. . (or  -G 
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G^jD  ) matrix  with  the  appropriate  dressing.  Similarly  for  the  ponderous 
matrices  we  have 
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(172) 


The  various  quantities  required,  in  addition  to  the  coefficients  in  the  original 
form  of  the  perturbation  equation  of  motion,  to  transform  to  another  form  of  the 
equation  are  stmmarised  in  Table  13.  The  transformations  between  the  various 
forms  derived  using  body-fixed  axes,  require,  when  there  is  a change  between 
encastrd  and  free-free  modes,  some  modal  information  in  the  form  of  the  matrices 
Rq,  . In  addition  certain  information  on  the  datum-motion  forces  is  required 
when  transforming  to  or  from  the  constant-velocity  axes  form.  For  example 
is  required  with  replaced  in  turn  by  Ae^,  ^pf»  ^sf 


(5  gdmAj^  ) . 


The  sum  of  the  forces  on  any  particle  during  the  datum  motion 


will  be  zero*  a^d  so,  for  example 
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One  can  therefore  assume,  in  transforming  to  or  from  the  constant-velocity  axes 
form,  that  all  the  individual  datum  motion  forces  are  zero.  The  result  will  be 
an  equation  which  has  the  right  solution  though  the  individual  terms  are  in 
general  wrong,  being  different  from  those  obtained  when  the  equation  is  derived 
directly.  Thus,  for  example,  for  a rigid  aircraft  (ie  n > 0)  in  a vacuum,  since 


from  the  equilibrium  equation  (Table  1 or  2)  we  then  have 
have  (of  Tables  3 and  4) 


-Ax, 


pf  ’ 


we 


and 


* This,  of  course,  is  more  than  is  said  by  the  datum  motion  equations  for  it  is 
inconceivable  that  the  degrees  of  freedom  should  be  as  numerous  as  the  particles 
of  the  aircraft. 
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in  the  constant-velocity  form  and  body-fixed  form  (displacement  body  freedoms) 
of  the  perturbation  equations  respectively.  However,  if  we  transform  from  one  to 
the  other,  using  equation  (161),  and  in  the  transformation  assume  the  datum  motion 
forces  are  zero,  then  we  find  that  these  gravitational  and  propulsive  matrices 
are  unchanged. 

7.4  With  certain  other  published  equations 

A number  of  authors  (Milne^,  Etkin^,  Taylor',  and  others)  have  suggested 
the  use  of  mean-body  axes  in  the  derivation  of  the  equations  of  motion  of  a 
deformable  aircraft.  In  our  derivations  of  sections  6.2.3  and  6.2.4  a frame  of 
reference,  occupying  the  location  the  body-fixed  axes  would  have  if  there  were  no 
deformation  (ie  if  the  first  n generalised  coordinates  were  zero) , can  be 
thought  of  as  mean-body  axes  provided  the  modal  matrix  R satisfies  the 
conditions*  for  mean-body  axes,  viz: 


^5mR  = 0 

[6mAx^R  ■ 0 


(177) 


The  equation  given  in  detail  in  Table  7 does  not  agree,  for  example,  as  regards 
individual  terms  with  that  obtained  by  Etkin^.  Like  us  he  used  Lagrange's 
equation  for  a non-inertial  frame  in  conjunction  with  body  freedom  equations 
obtained  by  the  principles  of  momentum.  However  he  naturally  considered  the 
components  of  the  momentum  referred  to  the  mean-body  axes  rather  than  the  body- 
fixed  axes.  It  is  easily  seen  that  the  mean-body  axes  become  coincident  with 
the  body-fixed  axes  if  they  are  given  translations 
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and  rotations 


(179) 


* A way,  but  not  the  only  way,  of  satisfying  (177)  is  to  take  as  modes  a set  of 
the  free-free  normal  mdes  in  the  absence  of  gravity. 
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(of  equations  (H5),  (2)  and  (M)).  If  we  denote  by  superscript  (m)  quantities 
referred  to  the  mean-body  axes*  we  can  write,  using  (5), 


* X - A= 


(180) 


L - (AxfRo 


(181) 


where  X,  L etc  can  be  written  in  terms  of  the  q^  using  equations  (93)  , (94) 
and  (112).  Thus  if  we  add** 


Ay  -P 
^gf  q 


0 0 


0 0 


AXgfRo  ^ ^LgfPq  0 0 


(182) 


to  [G. .] , and  apply  the  conditions  (177);  and  similarly  with  the  matrices  [P..] 

-j. 

and  [-Q. .] , we  should,  and  indeed  do,  get  the  body-freedom  equations  given  by 

. 7 , 

Etkin  . He,  incidentally,  follows  the  fairly  common  practice  of  including  the 
propulsive  terms  in  the  aerodynamic  matrix;  and  he  also  chooses  his  modes  to 
satisfy  the  normal  mode  condition,  that  ^6mR^R  is  diagonal,  in  addition  to  (177) 

It  also  can  be  shown  that  our  deformation  freedom  equations  of  Table  7, 
subject  to  the  mean-body  axes  condition  (equation  (177))  become  identical  with 
Etkins  , in  their  non-concise  form^,  if  a term  (of  equations  (161)  and  (163)) 


* Thus  t typical  moment  components  about  the  origin  of  the 

mean-body  axes  (hence  the  subscript)  and  referred  to  those  axes  (hence  the 
superscript) . 

**  In  fact  “ 0 ; and  there  also  is  no  change  to  lE^jl  since  the  structure 

cannot  exert  any  overall  force  on  itself. 

t In  the  'non-concise'  form  they  had  before  division  by  an  appropriate  mass  or 
moment  of  inertia. 
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(183) 


is  added  to  [^£j1  i-n  the  1 1 submatrix  position,  and  similar  additions  are  made  to 
[^ijl  • snd  • Some  of  these  additional  terms  are  of  course  zero  as  a 

consequence  of  the  characteristics  of  the  gravitational  and  structural  forces 
(sections  5.3  and  5.5)  and  of  the  conditions  (177). 

Sunming  up,  therefore,  we  see  that  we  can  transform,  from  our  perturbation 
motion  equation  obtained  using  body-fixed  axes,  free-free  modes,  and  velocity 
body  freedoms  (Table  7 and  section  6.2.3)  with  deformation  modes  satisfying  the 
conditions  (177),  to  equations  obtained  using  mean-body  axes  and  similar  free- 
doms by  adding 


- Ia  R - (Ia  a V 0 0 


(184) 
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to  and  making  similar  additions  to  [P..],  [E.  .]  and  [-Q..].  Since,  as 

ij  ij  ij 

noted  in  section  7.3,  the  sum  of  the  forces  on  any  particle  during  the  datum 
motior  will  be  zero,  it  follows  that  the  sum  of  these  additions  will  be  zero. 

The  resultant  equations  are  given  in  detail  in  Table  14.  Note  that  the  aero- 
d)niamic  overall  force  and  moment  about  the  origin  of  the  mean-body  axes,  referred 
to  the  mean-body  axes,  are  (<;/ equations  (155),  (156),  (180)  and  (181)) 
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Though  considering  the  use  of  various  frames  of  reference  - mean-body  axes, 
body-fixed  axes,  principal  axes  - Milne  (Ref  6,  Part  I,  Appendix  I)  states  the 
perturbation  equations  of  motion.  a concentrated  form  appropriate  to  any  type 
of  body  axes.  Thus  his  equations,  when  the  perturbations  are  small,  with  the  body 
axes  chosen  to  be  coincident  with  the  principal  axes  of  inertia  through  the  air- 
craft centre  of  gravity  during  the  datum  motion  (ie  equations  (52)  and  (53)  being 
satisfied),  are  identical  with  ours  for  body-fixed  axes,  free-free  modes  and 
velocity  body  freedoms  (Table  7 and  section  6.2.3)  when  the  matrices  (184)  etc 
are  added  to  our  constituent  matrices  of  Table  7. 

For  a rigid  aircraft  it  is  usual  to  reduce  the  perturbation  equations  of 
2 

motion  to  what  Hopkin  called  'a  concise  form'  by  dividing  each  force  equation  by 
the  mass  and  each  moment  equation  by  the  appropriate  moment  of  inertia.  Etkin^ 
and  others  (Ref  8 for  example)  have  similarly  set  up  a 'concise'  form  for  the 
equations  of  motion  of  the  deformable  aircraft.  Thus  their  final  form  ’ is 
that  given  when  the  equation  of  Table  14  is  premultiplied  by 


^6mR^R  0 0 

0 ml  0 

0 0 1 

ij 


(187) 


Finally  we  would  just  mention  that  various  other  authors  (Refs  9 to  13  et  al) 
have  derived  or  quoted  the  perturbation  equations  of  motion  for  a deformable  air- 
craft. For  the  datum  motion  that  we  have  considered  it  is  believed  that  their 
forms  are  all  included  in  the  cases  we  have  already  considered  though  the  nota- 
tion and/or  assumptions  used  are  not  always  clear. 

8 CONCLUDING  REMARKS  ♦ 

We  have  been  able  to  develop  in  detail  various  forms  of  the  perturbation 
equations  of  motion  of  a deformable  aircraft  and  present  them  in  the  attached 
tables  (Tables  3 to  7 and  14).  These  relate  to  small  perturbations  from  a datum 
motion  which  is  straight  flight  with  constant  linear  and  zero  angular  velocity. 

The  relationships  between  the  various  forms  were  established  and  stated,  and  the 
requirements  for  transformation  from  one  form  to  another  listed  (Table  13).  A 
particular  feature  of  the  work  has  been  the  care  with  which  account  has  been 
taken  of  the  various  forces  (structural  etc)  which  were  already  present  in  the 
unperturbed  state.  To  make  use  of  the  equations  the  reader  will  of  course  have 
to  call  in  the  help  of  the  aerodynamic is t and  the  elastician  (or  make  use  of 


i 


J 


077 


64 


the  results  of  a ground  resonance  test  as  in  Appendix  C)  and  possibly  others 
as  well. 

There  is  room  for  further  development:  different  unperturbed  conditions, 
better  representation  of  the  propulsive  forces,  more  realistic  model  of 
structural  damping,  more  accurate  representation  of  the  hereditary  nature  of  the 
aerodynamic  forces,  etc  - but  the  need  for  such  refinements  has  yet  to  be  con- 
vincingly demonstrated. 
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Appendix  A 

STRUCTURAL  FORCES  FOR  AN  ISOTROPIC  ELASTIC  BODY 


The  meaning  of  the  term  'structural  force'  on  a particle  or  element  is 
perhaps  better  understood  when  it  is  defined  in  terms  of  stresses  for  an  iso- 
tropic elastic  body  with  no  structural  damping.  Then  (of  eg  Ref  3)  the  structural 
force  is  given  by,  for  the  unperturbed  state 
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(A-1) 


where  a^  etc  are  the  stress  components  (a^  ■ etc) 

J means  the  integral  over  the  free  surface  of  the  element 
surface 

£>m,n  are  the  direction  cosines  of  the  outward  drawn  normal  to  the 
surface. 

The  first  term  is  the  structural  force  there  would  be  on  an  entirely  internal 

element  of  the  body,  and  the  second  is  the  reduction  consequent  from  the  fact  that 

there  is  no  structure  to  exert  a force  on  the  free  surface  of  an  element.  It  is 
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immediately  obvious,  using  Green's  formula  relating  surface  and  volume  integrals, 
that  the  sum  of  the  above  structural  forces  over  the  whole  body  is  zero  and  so 
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We  also  have  from  Green's  formula,  with  S enclosing  the  volume  V 
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Using  this  and  similar  relationships  derived  in  the  same  way  we  find  therefore 

I that 
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we  get  from  the  first  order  terms  the  relationship 
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The  latter  equation  is  of  course  equivalent  to  (A-7) . Now  the  general  expres- 
sions for  the  overall  forces  and  moments,  referred  to  the  constant-velocity 
axes,  (equations  (74)  and  (75))  give,  using  (A-8)  or  (46),  structural  contribu- 
tions due  to  the  perturbations; 
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Consequently  from  (A-2),  (A-9) , (A-13)  and  (A-14)‘all  these  contributions  are 
zero. 


The  above  expressions  (A- 15)  and  (A- 16)  are  of  course  independent  of  the 
assumption  of  an  isotropic  elastic  body;  and  so  if  we  assert  that  a structure 
cannot  exert  any  overall  translational  force  or  moment  on  itself  we  immediately 
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find  that  the  relationships  (A-2) , (A-7) , (A-9) , (A-13)  and  (A-14)  are  true  in 
general . 

So  far  we  have  taken  no  account  of  the  relationships  between  stress  and 
strain  and  displacement.  The  strain  components  are  simply  related  to  the  first 
derivatives  of  the  displacements;  and  for  our  'isotropic  elastic  body'  the  stress 
and  strain  components  are  linearly  related  by  Hooke's  Law.  Consequently  it  can  be 
shown  that  the  strain  energy  of  a deformed  elastic  body  is  a quadratic  function 
of  the  first  derivatives  of  the  displacements  (of  Ref  14).  Thus,  in  our  case, 
the  strain  energy  E is  a quadratic  function  of  the  generalised  coordinates  q.  : 
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where,  without  any  loss  in  generality,  lE^j]  is  written  as  a symmetric  matrix. 

From  the  fact  of  the  existence  of  such  a potential  function  it  immediately 
follows  that  the  matrix  of  the  generalised  structural  force  coefficients,  in  the 
perturbation  equation  of  motion,  is,  as  we  have  already  indicated  by  our  notation, 
[E^j]  which  is  symnetric.  Thus  from  the  form  already  deduced  for  this  matrix 
(Table  3)  we  see  that 
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PERTURBATION  EQUATICTJS  OF  MOTION  NEAR  THE  GROUND 

For  an  aircraft  in  flight  close  to  the  ground  the  boundary  condition  to  be 
satisfied  by  the  air  flow  at  the  surface  of  the  ground  has  an  important  effect  on 
the  air  forces  applied  to  the  aircraft.  Consequently  the  forms  taken,  in  the 
main  part  of  this  Report,  for  the  local  aerodynamic  force  vector  at  the  surface 
of  the  aircraft  (equations  (37),  (38),  (101),  (110)  and  (118))  are  no  longer  an 
adequate  representation  of  the  truth.  It  is  necessary  then  to  write  the  local 
aerod3mamic  force  vector  in  the  appropriate  one  of  the  forms  we  have  taken  for 
the  typical  local  force  vector  (of  for  example  equations  (120)  and  (121)  in 
section  A.l).  The  matrix  of  generalised  aerodynamic  force  coefficients,  for  the 
constant-velocity  axes  derivation,  is  consequently  given  by  equation  (78)  with 
the  dressing  omitted.  Similarly  the  other  forms  are  obtained  using  the 
relationships  (161),  (163),  (164)  and  (165)  - again  omitting  the  bars  -,  and 
put  in  terms  of  the  desired  local  and  overall  aerodynamic  load  coefficients  by 
the  use  of  equations  (126),  (130),  (144),  (145),  (146)  and  (147).  The  resultant 
aerodynamic  matrices  are  given  in  detail  in  Table  15.  For  the  various  body-fixed 
axes  cases  the  local  aerodynamic  force  vector  has  been  written  as  (of  equations 
(100),  (104)  and  (112)) 


n+6 
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^f 
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where 


1 • 


"n+6 
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l\ 


(B-2) 


and  the  coefficients  e ^ , e^  etc  are,  in  general,  differential  operators. 

Similar  expressions  hold,  and  have  been  used,  with  the  subscripts  (x,y,z,(^,0,i|)) 
replacing  (n+1 , . . . ,n+6) , for  the  overall  forces  and  moments.  Comparing  (37)  with 
(B-I),  and  so  on,  we  note  that 


'n+ 1 * ' ’ ®n+6 

®n+I  *■*  ®n+6 
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[fn+l  •_ 
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(B-3) 


etc,  as  the  distance  of  the  aircraft  from  the  ground  becomes  large. 

B. 1 Upholding  forces 

When  there  is  contact  between  the  aircraft  and  the  ground,  as,  for  example, 
during  the  take-off  run,  additional  forces  are  obviously  felt  by  the  aircraft. 
These  we  will  call  upholding,  or  support,  forces.  We  will  assume  the  normal 
earth-fixed  axes  Oq*o^O*0  been  defined  such  that  during  the  datum  motion 

the  aircraft  has  zero  velocity  components  in  the  directions  O^yQ  and  OqZq  . 

Its  velocity  u^  in  the  direction  O^Xq  will  be  constant,  in  accord  with  the 
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assumption  throughout  this  Report  that  any  datum  motion  considered  is  one  of 
zero  angular  velocity  and  constant  linear  velocity.  One  case  of  certain  particu- 
lar interest  (af  for  exaiiq)le  Appendix  C)  is  that  for  which  u^  is  zero.  We  also 
assume  the  earth  surface,  when  under  no  external  load,  is  given  by 


(B-4) 


and  that  the  earth  behaves  elastically*  such  that  when  displaced  to 
Zq  =•  z^^\xQ,yQ)  it  exerts  a force 


JO)  _ (0) 

z • z 

p - 


(B-5) 


at  (xQ,yQ,z''  '),  referred  to  the  normal  earth-fixed  axes. 

Transformation  from  the  normal  earth-fixed  axes  to  the  constant-velocity 
is  accomplished  by 

(i)  translations  (u^t,©,©)  in  the  direction  of  the  normal  earth-fixed 


(ii)  rotations  (^>^,0^,^^)  in  the  usual  Euler  order  about  the  carried  axes. 

Thus  we  are  assuming  that  the  reference  point  coincides  with  the  origin  0^  of 
the  normal  earth-fixed  axes  at  the  instance  t = 0 during  the  datum  motion. 
Writing 


S’ 


(B-6) 


we  find  that  a particle  on  the  outside  of  the  aircraft,  which  was  at  the  point 
(Xf,yf,Zf)  in  the  unperturbed  state,  is  in  contact  with  the  ground  at  an  instant 
during  small  perturbations  of  the  d.atum  motion  provided  that 


* If  desired  a more  general  representation  could  be  used  with  6'  ^ taken  to  be 
a differential  operator.  ^ 
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where  and  their  derivatives  are  evaluated  at  the  point  given  by 

(B-8) . 

Now  during  the  datim  motion  these  forces  must  be  constant  and  so  we  must 
have  either 


Uq  - 0 


(B-12) 


everywhere  . (B-13) 


With  either  of  these  conditions  satisfied  we  can  write,  at  the  point  given  by 
(B-8) , 
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Thus  (B-11)  becomes,  remembering  that  S is  orthogonal 
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Note  that,  from  (B-14), 
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(B-15) 
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and  so 


Similarly 
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Thus  the  functions  6 , z must  be  such  that 
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(B-22) 


In  addition,  if  the  condition  (B-13)  is  the  one  which  is  satisfied  during  the 
datum  motion,  then  these  two  functions  must  also  satisfy 
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(B-24) 


Expressing  the  perturbation  in  terms  of  the  generalised  coordinates  used  in 
the  constant-velocity  axes  development  (equation  (7)),  equation  (B-17)  becomes 
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(B-25) 


The  resultant  overall  forces  and  generalised  forces  can  be  obtained,  as  in  the 
other  cases,  from  the  general  expressions  (74),  (76),  (22)  and  (78).  We  will 
write  a generalised  upholding  (support)  force  as 


«i>f 


(B-26) 


and  the  perturbation  equation  of  motion  will  then  be 

{(A..l»^  . IS.. I . (G..1  . IP., I . 1E..1  - [Q..l| 


*n+6 


0 (B-27) 


where  the  are  obtained  as  above,  the  aerodynamic  matrix  iQ^jl  is  given  in 

Table  15  and  the  other  matrices  are  given  in  Table  3.  The  equation  of  equilibrium 
(Table  1)  can  similarly  be  modified. 
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Appendix  C 

PERTURBATION  EQUATIONS  OF  MOTION  IN  THE  GROUND  RESONANCE  TEST  SITUATION 


During  a ground  resonance  test  it  is  reasonable  to  assume  that  there  will 
be  no  aerodynamic  or  propulsive  forces,  and  the  constant-velocity  axes  are 
stationary  and  such  that  O^y^  is  coincident  with  the  normal  earth-fixed  axis 
O^yQ  . Thus  the  two  attitude  angles  of  the  undisturbed  aircraft  are  zero. 

There  will  of  course  normally  be  some  excitation  forces  but  as  throughout  the 
rest  of  this  Report  we  will  assume  them  to  be  zero.  It  will  be  a simple  matter 
for  the  reader  to  include  them  if  desired.  We  will  also  assume  that  at  each 
support  point  the  ground  profile  and  stiffness  are  both  stationary.  Hence  from 
(B-25)  the  upholding  (support)  force  vector*  is,  referred  to  the  constant- 
velocity  axes 
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where,  in  this  ground  resonance  test  situation 
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From  the  equilibrium  equation  (.of  Table  1)  we  immediately  find  that 
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and 


S *f^f 


(C-5) 


* We  are  continuing  to  assume  that  the  support  forces  are  purely  elastic.  For 
the  soft  suspension  desirable  for  a ground  resonance  test  it  may  be  necessary 
to  take  account  also  of  the  effect  of  the  mass  of  the  suspension  on  these 
support  forces. 
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and  so,  using  equations  (75),  (77)  and  (78)  along  with  (C-l),  the  upholding 
(support)  matrix  is  given  by 


0 -g(j:6mR^)AA^  - I (C-6) 


gAo  j;6mR  + lAe  fR  0 


where  we  have  made  use  of  the  argument:'- (C-3)  is  true  for  any  R and  in 
. T T 

particular  for  R = R A^  where  A^  is  an  arbitrary  constant  skew-synmetric 
matrix  of  the  type  of  equation  (6) , and  so 
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which  gives  since  {a  ...}  is  arbitrary 


■ (■<’■') 


The  gravitational  and  structural  matrices  are  given  in  Table  3,  and  so 
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The  equation  of  motion  for  small  perturbations  is  therefore 


+ [S..I  + [G..] 
ij  ij 


(C-10) 


where  the  inertia  matrix  is  given  in  Table  3 and  the  other  terms  come  from  the 
above  equation  (C-9) . It  is  normal  practice  in  ground  resonance  tests  to  make  the 
support  stiffnesses  (the  5^)  as  small  as  possible  so  that  the  first  term  in 
(C-9)  can  generally  be  ignored.  We  can  also  replace  the  terms  involving  the 
datum  values  of  the  structural  force  vectors  by  terms  involving  the  coefficients 
in  the  expressions  for  the  perturbations  thereto.  This  is  done  through  the 
relationship 


(C-ll) 


which  is  a consequence  of  the  fact  that  the  structure  cannot  exert  any  overall 
moment  on  itself  (of  the  particular  demonstration  in  Appendix  A) . It  also  should 
be  noted  that  for  most  structures  (of  Appendix  A)  the  submatrix  in  (C-9) 
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can  be  taken  to  be  symmetric. 


The  main  objective  of  a ground  resonance  test  frequently  is  to  gather 
information  about  the  structural  forces.  If  the  function  R represents  the 
normal  modes  measured  in  such  a test  which  are  not  predominantly  body  freedom 
modes • and  the  suspension  is  sufficiently  soft  then  we  see  that  to  a good 
approximation 

^6mR''^R  ■ is  diagonal  (C-12) 


^6mR  “ 0 


(C-13) 
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is  diagonal . 


(C-14) 


(C-15) 


Thus  knowing  the  frequencies  of  these  normal  modes  we  can  immediately  find  all 
the  elements  in  the  second  term  in  (C-9) , or  in  the  structural  matrix  of  Table  3, 
for  these  particular  degrees  of  freedom  and  this  particular  datum  motion.  That 
is  we  have 


diagio) 
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where  is  the  circular  natural  frequency  of  the  ith  mode.  We  would  emphasise 

that  these  values  of  the  structural  coefficients  strictly  apply  only  to  the 

ground  resonance  test  situation.  For  other  stress  states  during  the  datum 

motion  it  is  obvious  from  (C-11)  or  (A-13)  that  fe  , ...  e "I  will  be  different 

sn  1 


even  if  the  modal  matrix  R is  the  same.  Nevertheless,  being  devoid  of  hope  of 
doing  anything  better  without  excessive  effort,  it  will  often  be  assumed  that 
these  coefficients,  e . etc,  in  the  components  of  the  local  structural  forces, 
are  independent  of  the  chosen  unperturbed  condition. 

To  obtain  the  structural  matrix  in  the  perturbation  equations  of  motion 
using  body-fixed  axes  is  a little  more  difficult  (of  Tables  4,  5,  6 and  7).  Wt 
have  from  (131)  and  (46): 
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which,  for  normal  modes,  can  be  evaluated  from  (C-II),  (C-16)  and  (C-17).  The 
other  terms  for  example,  Table  4)  comprise  one,  * which  similarly 

can  be  easily  evaluated,  and  one  which  sets  a problem.  To  be 

able  to  evaluate  the  latter  it  appears  that  one  will  also  require  measurements  of 
the  unperturbed  values  of  the  upholding  (support)  forces.  When  using  constant- 
velocity  axes  it  was  remarked  above  that  one  may  often  have  to  assume  that 


sn 


, and  hence  (c/,  for  example,  (C-ll))  the  weighted  sum  J^Ae^^R  are 
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independent  of  the  chosen  unperturbed  condition.  Similarly  with  body-fixed 
axes  one  may  have  to  assume  that 
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sn 
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independent  of  the  unperturbed  state. 

In  section  7.4  we  considered  also  the  use  of  mean-body  axes.  It  is 
however  rather  more  difficult  to  relate  the  structural  coefficients  there 
derived  (.of  Table  14)  to  the  results  of  a ground  resonance  test,  for  the  normal 
modes  obtained  in  a ground  resonance  test  do  not  satisfy  the  two  conditions 
(equation  (177))  for  mean-body  axes  (c/ equations  (C-H)  and  (C-14)).  We  will 
not  therefore  pursue  the  connection. 


This  consideration  does  however  highlight  the  fact  there  is  in  general 
some  inertia  coupling  between  the  normal  modes  measured  in  a ground  resonance 
test  and  the  rigid  body  rotational  freedoms.  How  large  this  is  it  is  difficult 
to  tell.  No  doubt,  for  an  almost  two-dimensional  structure,  such  as  an  air- 
craft, it  could  certainly  be  made  small  by  the  use  of  a continuous*  soft  support  - 
from  (C-II)  and  (C-14)  it  clearly  will  be  small  if  the  structural  forces  in  the 
unperturbed  state  are  small.  In  practice  the  cross  inertia  couplings  between 
the  measured  modes  and  the  rotational  body  freedoms  have  rarely  been  calculated. 
Moreover,  if  they  had  been,  one  would  not  know  how  much  was  the  'true'  value, 
how  much  was  due  to  experimental  error  (it  is  not  unusual  to  find  large  inertia 
couplings  between  two  measured  'normal'  modes),  and  how  much  was  due  to 

* Continuous  in  the  sense  that  every  point  on  the  underside  of  the  aircraft  is 
supported. 
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inadequacies  in  our  mathematical  model  of  the  suspension  (massless,  elastic, 
and  infinitely  soft).  All  this,  along  with  considerations  of  structural  damping 
and  friction,  suggests  that  measured  normal  modes  should  be  used  with  caution, 
and  that  one  should  regard  ground  resonance  tests  primarily  as  a fair  check  of 
theoretical  predictions  rather  than  a source  of  data. 
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Table  2 

DATUM  MOTION  EQUATICBi  USING  BODY-FIXED  AXES 

[-  (Vj  ♦ cepj  * (Pj),  * (ipj  - 0 


Constituent  matrices 


Aerodynamic 

Gravitational 

Propul a ive 

Structural 

- i(p.),i 

- '■JVPfl  1 I 


K - W 


Overall  forces 


See  Table  I 


Note;  For  the  purpose  of  this  table  equations  (29),  (30)  and  (31)  have  been 
combined  into  one  matrix  equation. 
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Table  3 

PERTURBATIOM  MOTION  EQUATION  USING  CONSTANT-VELOCITY  AVF.R 


Equation 
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Constituent  matrices 

(i)  Inertia  (ie  ponderous  inertia) 


lAi.l 


^6inR  R ^6mR  -^6inR^Ax£ 
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(ii)  Gravitational 
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Table  4 

PERTURBATION  MOTION  EQUATION  USING  BODY-FIXED  AXES  IN  TERMS  OF  ENCASTRE  MODES 

AND  DISPLACEMENT  BODY  FREEDOMS 
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(iii)  Propulsive 
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Table  4 (concluded) 


Constituent  matrices  (concluded) 
(iv)  Structural 
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Overall  aerodynamic  force  coefficients 

See  Table  9 
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Table  5 


PERTURBATION  MOTION  EQUATION  USING  BODY-FIXED  AXES  IN  TERMS  OF  ENCASTRE  MODES 

AND  VELOCITY  BODY  FREEDOMS 


Equation 


* * rtiji  ♦ IC..1  ♦ (f..l  ^ 1E..1  - IQ, 


iQjji  * [G|.iQ-'}  r , 


»a> 


90 


Table  5 (continued) 


Constituent  matrices  (continued) 
(iii)  Propulsive 
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Table  6 


PERTURBATIOM  MOTION  EQUATION  USING  BODY-FIXED  AXES  IN  TP.BMr  m?  i?pbb-free  MODES 

AND  DISPLACEMENT  BODY  FREEDOMS 


Equation 


{K..)d2  . By)  * [Pyl  * ffy,  - By,} 
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Constituent  matrices 

(i)  Ponderous  inertia 
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(ii)  Gravitational 
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(iii)  Propulsive 
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0 0 
0 0 


^lAcpfR  + (lAcpfAx^jp^  - AxpfRg 


0 0 


Table  6 (concluded) 


Constituent  matrices  (concluded) 
(iv)  Structural 


Overall  aerodynamic  force  coefficients 


See  Tables  9 and  1 1 
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Table  7 


PERTURBATICTt  MOTIOM  EQUATICT}  USING  BODY-FIXED  AXES  IN  TEBMS  OF  FREE-FUgw  MODES 

AND  VELOCITY  BODY  FREEDOMS  ~ ^ 


Equation 


{k..)d2  . [I.. Id  ♦ [5j.l  * [6j.|  * CPj.1  « (S..1  - (q^.l  * (g,.|D-'}  h 1 . o 


Constituent  matrices 


(i)  Ponderous 


IdmR^R  0 0 
J]6mR  0 0 


^ 6mAx^R  ^ ^ 


0 -^dmR'^Axj 


0 0 


0 0 


0 0 


(ii)  Gravitational 


%P-^)  - 


*Xgf“o 


95 


077 


Table  7 (continued) 


Constituent  tnatrices  (continued) 

(ii)  Gravitational  (concluded) 


le,.) 


0 

0 

0 


An 


-Ax, 


gf 


(iii)  Propulsive 


if..i 


0 


^AcpfR  + (^^pf^Xf)^q  ^XpfRg 


0 0 
0 0 
0 0 


(iv)  Structural 

-Ir 


g,.i  - 


Trj 


si 


sn 


’si 


^qj^^esfR  + (^^egfAx^jPqj 


0 0 


0 0 


0 0 


J 


r*) 


f 
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T 


K 


; j 

i : 


il- . 


[£ 
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Table  9 


OVERALL  FORCE  COEFFICIBMTS  - BODY-FIXED  AXES  ••  ENCASTRE  MODES 
AND  DISPLACEMENT  BODY  FREEDOMS 


Gravitational 


^ * 

gl 


g> 


2 1 
L8‘ 


gl 


M 


gl 


|N  I 

L gl 


Propulsive 


^ « 

pi 


pi 


L pi 


pi 


M 


pi 


N 


pi 


g.n+6 


[0  0 Ajcgfl 


g,n+6 


^XgfRo  ■ I^tdR  0 oj 


p,n+6 


p,n+6 


Aerodynamic  - well  away  from  ground 


X.  ...  X 

1 n 

- I 

e,  ...  e 

I n 

Y 

f 

1 

1 

A 

2 

Q ......... 

Lfi 

[fl  J 

Table  9 (concluded) 


where  X-  - . 

H * **=*'• 

A 2 

NB:  X*  , for  exantple,  is  not  the  usual  X;  derivative.  All  the  above  terms 
such  as  X*«  X*  are  in  general  differential  operators,  and  thus  there  will  be 

• “ • A A A 

contributions  to  the  X?  derivative  from  X*.  X*  and  X*  . 
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Table  10 

OVERALL  AERODYNAMIC  FORCE  COEFFICIENTS  FOR  VELOCITY  BODY  FREEDOMS  - 
BODY-FIRED  AXES  - WELL  AWAY  FR(»1  GROUND 
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Table  12 

OVERALL  AERODYMAMIC  FORCE  COEFFICIENTS  - BODY-FIXED  AXES  - PREE-pggg 

VELOCITY  BODY  FREEDOMS  - AWAY  FROM  GROOin) 


X|  • • • X 
1 n 


Zfe,  « 


L,  ...  t 
1 n 


1 


...  e. 


“ lAe^R  + AxjRq  - ^jAejA*jjp^ 


For  other  aerodynamic  coefficients  see  Table  10  since  X ■ X . X - X etc. 

u u p p 
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^ — •Wf 
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Table  13 


TRANSFOBMATION  BETWEEN  DIFFERENT  FORMS  OF  THE 


PERTDRBATIOl  EQUATION  OF  MOTION 


To  transform  from  one  form  of  the  perturbation  equation  to  another  one 
requires,  in  addition  to  the  coefficients  in  the  original  equation,  the  following 
quantities  (of  section  4.3): 


Body-fixed 

To 

(or  from) 

Constant- 

Encastr£ 

Free-free 

From 
(or  to) 

velocity 

Disp. 

Velocity 

Disp. 

Velocity 

Dressing 

Form 

Dressing 

A 

w 

Constant- 

velocity 

Body- 

fixed 

Encastr^ 

displaceaient 

A 

V % 

Xj,  etc 

I*ef®» 

1 

Body- 

fixed 

EncastrC 

velocity 

A 

*^0’  ***f 

Xj,  tj,  etc 

Auf 

Body- 

fixed 

Free-free 

displacement 

etc 

"q 

Auf 

i 

Body- 

fixed 

Free-free 

velocity 

W 

*0’  ^q’  *“f 

X£t  etc 

Z^ej^Xf 

’^0'  ^ 
Au£ 

^)*  '■q 

A 

"f 

Notes:  (i) 

(ii) 

(iii) 

(iv) 


Au^  should  always  be  known. 

Need  information  about  deformation  modes  going  from  or  to  to  deduce 
^q  • 

Entries  in  the  table  of  a typical  force  vector  indicate  the  complete 
set  of  these  (aerodynamic,  propulsive  etc)  are  required. 

In  transformation  from  or  to  the  constant-velocity  form  the  datum- 
motion  forces  could  be  considered  zero.  The  result  would  be  an 
equation  with  the  wrong  constituent  coefficients  but  the  right 
solution. 


106 


I : 

M 

I : 
r I 

LI 


1 


r I 


f-  : 
i 

t 

I '■ 

• j 


I i 


...jle  14 


PEBTURBATIOM  MOTION  EQUATIOM  USIHG  MEAH-BODY  AXES  IN  TERMS  OF 
FREE-FBEE  MODES  AND  VELOCITY  BODY  FREEDOMS 


Equation 


* OjjlD  ♦ IV.. 1 ♦ ff..I  ♦ Byl  - (5.. I * liSfjiD-'} 


Constituent  matrices 


(i)  Ponderous 


1X..1 


diagUfinR  R 0 0 


[J..] 

ij 


liagjj; 

liagjo 


diagfO  ml  I 


4 


rv,,i 


0 0 
0 0 


0 0 


0 

■®^f 

0 


(ii)  Gravitational 


I5*,l 


0 0 

0 0 -A: 

0 0 


«gf 

0 


(iii)  Propulsive 


Sj.i 


0 

0 


L’oHJ 


0 0 
0 0 


jAepfR  + (XAxfAepJPq  0 oj 


Table  14  (concluded) 


Constituent  matrices  (concluded) 
(iv)  Structural 


lEi-]  - 


'si 

A 

^8l 

A 

L^sl 


0 

0 


sn 


0 0 
0 0 


(v)  Aerodynamic 


...  s 

1 n 

» * 

®n+l  *'•  ®n+3 

6 i • • • 6 ^ 

n+4  n+6 

f 

f 

f 

1 

n+1 

n+4 

g 

i 

Q 

Lfi  j 

lrn+1  ll 

_n+4  _ 

"(c) 

• • • A 

1 n 

v(c) 

^ ^ 

XXX 

U V w 

Y 

^ ^ ^ 

XXX 
p q r 

A 

Y 

1 

u 

z 

P 

7 



L.1  J 

_ u 

L p J 

(c)  (c) 

L,  • • • hjj 

^ ^ 

L L L 

U V w 

fi  

^ A 

L L L 

p q r 

M 





u 

A 

N 

p 

N 

Ll  J 

L u J 

L p J 

NB:  In  the  above  equations  the  modal  matrix  R has  been  chosen  such  that 

^6mR  - 0 

^dmAxjR  ■ 0 

These  equations  of  motion  have  not  been  derived  directly  in  the  present  paper  but 
in  section  7.4  we  have  merely  shown  that  they  are  equivalent,  for  a modal  matrix 
satisfying  the  above  conditions,  to  the  equations  of  Table  7. 


c 

(i)  Fundamental  terms 
degree  of  freedom 


semi-rigid 


body  freedom 

perturbation 
deformation 
circular  frequency 
natural  frequency 

generalised  coordinates 
generalised  force 
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GLOSSARY  OF  TERMS 


a system  is  said  to  have  n degrees  of  freedom 
if  its  position  and  shape  can  be  uniquely  specified 
by  n parameters,  but  not  by  less  than  n 

a body  is  said  to  be  semi-rigid  if  it  has  a finite  [ 

number  of  degrees  of  freedom  in  addition  to  any 
body  freedoms 

a body  is  said  to  have  n body  freedoms  (n  < 6) 
if  its  position,  with  a given  shape,  can  be 
uniquely  specified  by  n parameters  but  not  by 
less  than  n 

a disturbance  from  a datum  state 
a perturbation  in  shape 

the  frequency  of  oscillation  multiplied  by  2ir 

a constant  frequency  of  free  oscillation  of  a 
system 

the  set  of  parameters  used  to  specify  the  position 

and  shape  of  a system  j 

the  coefficient  of  the  increment  of  a generalised 
coordinate  in  the  expression  for  the  virtual  work 
of  a system 


(ii)  Types  of  force 
aerodynamic* 

propulsive* 

gravitational 

structural 

upholding  (or  support) 

ponderous 

inertia 

damping 


exerted  by  the  air  on  the  external  surface  of  a 
body 

produced  by  the  operation  of  a propulsive  unit 
such  as  a jet  engine 

resulting  from  the  gravitational  attraction 
between  each  particle  and  the  earth 

resulting  from  the  stresses  in  a body 

resulting  from  contact  between  a body  and  the 
ground 

a reversed  effective  force  consequent  to  a system 
having  mass 

a force  which  is  proportional  to  the  second 
derivative,  with  respect  to  time,  of  a generalised 
coordinate 

a force  which  is  proportional  to  the  first 
derivative,  with  respect  to  time,  of  a generalised 
coordinate 


077 


* The  distinction  between  aerodynamic  and  propulsive  is  not  entirely  clear-cut. 
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(ii)  (concluded) 


stiffness 


i 


GLOSSARY  OP  TERMS  (continued) 


s force  idiich  is  proportional  to  a generalised 
coordinate 


(iii)  Frames  of  reference  (all  right-handed  orthogonal*  cartesian) 


normal  earth-fixed  axes 


constant-velocity  axes 


body-fixed  axes 


mean-body  axes 


no-deformation-body- 
fixed  axes 


datum-attitude  earth 


(iv)  Orientation 

angle  of  inclination 


angle  of  bank 


axes  fixed  relative  to  the  earth  with  the  z-axis 
vertically  downwards 

axes  having  constant  linear  and  angular  velocity 
relative  to  an  inertial  frame.  (The  ones  used  in 
this  Report  have  zero  angular  velocity.) 

(For  the  datum  motion  considered  in  this  Report 
these  are  always  taken  to  be  coincident  with  the 
datum-attitude  earth  axes  - see  below.) 

axes  whose  origin  and  orientation  are  fixed  in  a 
small  material  portion  of  the  body  which  is  such 
that  the  axes  always  remain  mutually  perpendicular. 
(The  ones  used  are  such  that  the  origin  coincides 
with  the  centre  of  gravity  of  the  body,  and  the 
axes  are  parallel  to  the  principal  axes  of 
inertia,  during  the  datum  motion.) 

axes,  with  origin  at  the  centre  of  gravity,  of 
the  body,  orientated  so  that  the  kinetic  energy 
relative  to  the  axes  is  a minimum. 

an  axes  system,  arbitrary  except  inrso-far  as  it 
is  of  the  same  order  of  nearness  to  the  body- 
fixed  axes  and  the  datum-motion-body-fixed  axes, 
with  respect  to  which  it  is  assumed,  in  the 
constant-velocity  axes  derivation,  that  the  posi- 
tion of  any  particle  of  the  body  is  given  as  a 
linear  function  of  deformational  generalised 
coordinates. 

axes  with  which  the  body-fixed  axes  coincide 
during  the  datum  motion.  (These  are  not  in 
general  body-fixed  axes.  They  are  earth-directed 
not  earth-fixed.) 


angle  between  the  x-axis  of  the  body-fixed  axes 
and  the  horizontal  plane 

angle  between  the  z-axis,  of  the  body-fixed  axes, 
and  the  vertical  plane  containing  the  x-axis  of 
the  same  frame 


* With  one  of  the  assumed  forms  of  deformation  (equation  (1))  the  body-fixed 
axes  only  remain  mutually  at  right  angles  for  small  perturbations  of  the  datum 
motion. 


(iv)  (concluded) 

nose-azimuth  angle 


Cv)  Deformation 

free-free  modes 


encastrd  modes 


GLOSSARY  OF  TERMS  (concluded) 


angle  between  the  projection  of  the  x-axis,  of 
the  body-fixed  axes,  on  the  horizontal  plane,  and 
the  x-axis  of  the  normal  earth-fixed  axes 


unconstrained  modes  of  deformation  - actually  a 
slight  constraint  has  been  imposed  in  that  the 
body-fixed  axes  are  kept  orthogonal  during  small 
deformations 

modes  of  deformation  constrained  to  have  zero 
values  and  zero  slopes  at  the  reference  point 


In  this  list  i,  j are  used  as  typical  dumay  subscripts. 


C 

ex 


I .1  . 
X y 


1 

z 


I 

n 


J. 

1 


J.  . 

IJ 

etc 

♦ 


etc 


skew-symmetric  matrices  involving  6,  i|>  , etc  (see  equation  (6)) 
ponderous  inertia  coefficient 

lower  triangular  matrix  involving  6,  ip  (see  equation  (56)) 

lower  triangular  matrix  involving  e,  f,  g,  x,  y,  z (see  equation 
(71)) 

differential  operator  d/dt 
strain  energy  when  no  perturbation 

generalised  structural  force 

structural  stiffness  coefficient 

generalised  gravitational  force 

gyros tatic  force 

gravitational  stiffness  coefficient 

gravitational  'build-up'  coefficient  (<?/  for  example  Table  5) 
unit  matrix 

principal  moments  of  inertia  of  undeforroed  aircraft 

= diag{I  I I } 

X y z 

a certain  coupling  force  between  the  rotational  body  freedoms  and 
the  deformational  freedoms  in  Lagrange's  equations  referred  to  a 
non-inertial  frame 

ponderous  damping  coefficient 

matrices  formed  from  the  elements  of  {(^  6 etc  (see  equation 
(57a)) 

aerodynamic  rolling  moment 
total  rolling  moment 

gravitational,  propulsive,  upholding  (support)  rolling  moments 

aerodynamic  rolling  moment  coefficients 


M 

M 


aerodynamic  pitching  moment 
total  pitching  moment 
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M ,M  ,M 
g P u 

M .M;,M  , 
X X 9 

M;,M  ,M„ 

♦ u’  p 


0 X y z 
c c c c 

0 *_y  S'-. 

m nr  m m 
°nVn^n 

°o*oyo^o 


Rj  »R2»^3 


LIST  OF  SYMBOLS  (continued) 

gravitational,  propulsive,  upholding  (support)  pitching  moments 

aerodynamic  pitching  moment  coefficients 

aerodynamic  yawing  moment 
total  yawing  moment 

gravitational,  propulsive,  upholding  (support)  yawing  moments 


aerodynamic  yawing  moment  coefficients 


position  of  reference  point  (particle) , origin  of  body-fixed  axes 

position  of  reference  point  when  perturbations  involve  no  deforma- 
tion, origin  of  no-deformation-body-fixed  axes  {of  section  3) 

body-fixed  axes 

constant-velocity  axes  (identical  with  datum-attitude  earth  axes) 
mean-body  axes 

no-deformation-body-fixed  axes 
normal  earth-fixed  axes 
generalised  propulsive  force 


matrix  of  modal  slopes  at  reference  point  (equation  (11)) 
propulsive  stiffness  coefficient 


generalised  aerodynamic  force 
total  generalised  force 

matrix  relating  angular  velocities  and  orientation  (of  equation  (36)), 
see  Ref  1 (Fart  II)  or  Ref  2) 

aerodynamic  coefficient 

modal  matrix  (equation  (2)) 

value  of  modal  matrix  at  reference  point 

certain  matrices,  see  equations  (59)  to  (61),  formed  from  the 
modal  matrix  elements 

axes  transformation  matrix  (of  equation  (5) , see  Ref  1 (Part  II) 
or  Ref  2) 

S where  necessary  to  designate  arguments  - in  this  case  6^, 
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-S. 

1 


ij 

T1.T2.T3 


V. . 

IJ 

w 

X 

X 

X .X  ,x 
g P u 


X^.Xi.X^. 

etc 

Y 

Y 

T .T  .Y 
g P u 


etc 

Z 

Z 

Z^.Z  ,Z 
g p’  u 


Z*.Z  .Z^ 
♦ u’  p 

etc 


LIST  OF  SYMBOLS  (continued) 

generalised  upholding  (support)  force 

upholding  (support)  stiffness  coefficient 

certain  matrices  whose  elements  are  particular  datum  motion 
particle  coordinates  (see  equations  (62)  to  (64)) 

centrifugal  potential  function  used  in  Lagrange's  equations 
referred  a non-inertial  frame 

ponderous  stiffness  coefficient 

kinetic  energy  relative  to  reference  frame 
overall  aerodynamic  force  resolute 
total  overall  force  resolute 

gravitational,  propulsive,  upholding  (support)  overall  force 
resolutes 

aerodynamic  force  resolute  coefficients 

overall  aerodynamic  force  resolute 
total  overall  force  resolute 

gravitational,  propulsive,  upholding  (support)  overall  force 
resolutes 

aerodynamic  force  resolute  coefficients 

overall  aerodynamic  force  resolute 
total  overall  force  resolute 

gravitational,  propulsive,  upholding  (support)  overall  force 
resolutes 

aerodynamic  force  resolute  coefficients 


i 


'.y  I 


a. 

1 


c. 

1 


e 


e. 

1 


x-component  of  modal  function  (equation  (2)) 
y-component  of  modal  function  (equation  (2)) 
z-component  of  modal  function  (equation  (2)) 
x-component  of  local  aerodynamic  force  vector 

coefficient  in  x-component  of  local  aerodynamic  force  vector  when 
near  the  ground 
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e ,e  etc 
g P 

e .>e  . etc 
gl  81 


etc 

g P 

£ .*f  . etc 
gl  81 


gg.Bp  etc 


LIST  OF  SYMBOLS  (continued) 

coefficient  in  x-conponent  of  local  aerodynamic  force  vector  when 
well  away  from  the  ground 

x-conq>onent8  of  gravitational,  propulsive,  etc  local  force  vectors 

coefficients  in  x-components  of  gravitational,  structural  etc 
local  force  vectors 

y-component  of  local  aerodynamic  force  vector 

coefficient  in  y-component  of  local  aerodynamic  foree  vector  when 
near  the  ground 

coefficient  in  y-conq)onent  of  local  aerodynamic  force  vector  when 
well  away  from  the  ground 

y-coiiq>onent8  of  gravitational,  propulsive  etc  local  force  vectors 

coefficients  in  y-conponents  of  gravitational,  structural  etc 
local  force  vectors 

acceleration  due  to  gravity 

z-component  of  local  aerodynamic  force  vector 

coefficient  in  z-component  of  local  aerodynamic  force  reactor  when 
near  the  ground 

coefficient  in  z-component  of  local  aerodynamic  force  vector  when 
well  away  from  the  ground 

z-components  of  gravitational,  propulsive  etc  local  force  vectors 

coefficients  in  z-components  of  gravitational,  structural  etc 
local  force  vectors 

first  column  of  (equation  (B-9)) 

second  column  of  (equation  (B-10)) 

direction  cosine  of  outward  drawn  normal  to  surface 
third  column  of  S^^  (equation  (39)) 

direction  cosine  of  outward  drawn  normal  to  surface 
mass  of  aircraft 

direction  cosine  of  outward  drawn  normal  to  surface 

number  of  deformational  degrees  of  freedom 

angular  velocity  resolute 

angular  velocity  resolute 

generalised  coordinate 

angular  velocity  resolute 

time 

linear  velocity  resolute 
particle  velocity  resolute 


j 

lJ 
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list  of  SYHBOLS  (continued) 


V linear  velocity  resolute 

V particle  velocity  resolute 
n 

w linear  velocity  resolute 

particle  velocity  resolute 

X particle  position  resolute 

Xj  resolute  of  reference  point  position  relative  to  its  position  in 

datum  motion 

X particle  position  resolute  relative  to  origin  of  constant-velocity 

axes 

X particle  position  resolute  relative  to  origin  of  mean-body  axes 

n 

X particle  position  resolute  relative  to  origin  of  no-deformation- 

body-fixed  axes 

y particle  position  resolute 

yj  resolute  of  reference  point  position  relative  to  its  position  in 

datum  motion 

y particle  position  resolute  relative  to  origin  of  constant-velocity 

axes 

y^  particle  position  resolute  relative  to  origin  of  mean-body  axes 

y particle  position  resolute  relative  to  origin  of  no-deformation- 

body-fixed  axes 

z particle  position  resolute 

Z|  resolute  of  reference  point  position  relative  to  its  position  in 

datum  motion 

z particle  position  resolute  relative  to  origin  of  constant-velocity 

axes 

particle  position  resolute  relative  to  origin  of  mean-body  axes 

z particle  position  resolute  relative  to  origin  of  no-deformation- 

body-fixed  axes 

ZpiXf.Yf.Zf)  “ Zp 

z^^^(x.,yn)  value  of  z-  at  the  surface  of  the  ground  when  under  no  external 
P ° ° load  ° 


angle  of  inclination 
angle  of  bank 

nose-azimuth  angle  (also  known  as  the  heading  or  heading  angle) 


LIST  OF  SYMBOLS  (continued) 


o etc 

xy 


■ (x  V ) 

earth  (or  support)  stiffness  coefficient 
mass  of  a particle 

orientation  angle  of  body-fixed  axes  relative  to  datum  motion 
stress  components 

orientation  angle  of  body-fixed  axes  relative  to  datum  motion 
orientation  angle  of  body-fixed  axes  relative  to  datum  motion 
circular  natural  frequency 

see  Table  IS 


Dressings 


(i)  Subscripts 

Absence  of  subscript,  in  appropriate  cases,  denotes  values  relative  to  or 
about  the  origin  of  the  body-fixed  axes. 

c denotes  value  relative  to  or  about  the  origin  of  constant- 

velocity  axes 

f denotes  values  during  datum  motion 

g denotes  gravitational 

m denotes  value  relative  to  or  about  the  origin  of  mean-body  axes 

n denotes  value  relative  to  or  about  the  origin  of  no-deformation- 

body-fixed  axes 

p denotes  propulsive,  or  (in  Appendices)  ground  profile 

s denotes  structural 

u denotes  upholding  (support) 

0 (nought)  denotes  value  relative  to  or  about  the  origin  of  normal  earth- 

fixed  axes 

(ii)  Superscripts 

Absence  of  a bracketed  superscript  to  a resolute  indicates  that  it  is  the 
value  of  the  resolute  along  the  body-fixed  axes. 

(c)  denotes  value  of  resolute  along  the  datum-attitude  earth  axes 

(ie  the  constant-velocity  axes) . This  is  a change  from  the  (f ) 
of  Ref  1 

(m)  denotes  value  of  resolute  along  the  mean-body  axes 
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